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ABSTRACT 


The solution to the steady state magnetohydrodynamic equations 
governing the supersonic expansion of the solar corona into interplane” 
tary space is obtained for various assumptions regarding the form in 
which proton thermal energy is carried away from the sun. 

The one- fluid, inviscid, formulation of the MHD equations is 
considered first assuming that thermal energy is carried away by con- 
•duction from a heat source located at the base of the corona. The 
inclusion in the £inalysis of the angular motion of the solar wind, leads 
to the existence of three critical points through which the numerical 
solutions must pass to extend from the sun's surface to large heliocentric 
distances. The results show that the amount of magnetic field energy 
converted into kinetic energy in the solar wind is only a small fraction 
of the total expansion energy flux and has little effect upon the final 
radial expansion velocity. 

The azimuthal velocity predicted by this model at 1 A.U. is 
tfl.l9 Km/sec., which is smaller that that indicated by experimental 
observations but in agreement with previous theoretical work in this 
field. 

The two-fluid formulation of the MHD equations is obtained next 
under the assumption that the protons become collisionless and thermally 
Beyond a gxven radius . This formulation is then applied to 
a two— region model of the solar wxnd xn which the flow in the inner 
region is described by the one- fluid equation?- and in the outer region 



by the two- fluid formulation. It is shown that the effect of the proton 
thermal anisotropy upon the angular motion of the solar wind is small and 
cannot increase the predicted azimuthal velocities at 1 A.U. to values in 
better agreement with observations. Since a modified CGL theory is used 
in the two- fluid formulation of the magnetohydrodynamic equations, the 
model provides, in addition, microscopic information about the protons in 
the form of velocity distribution function plots at various selected 
heliocentric distances. 

The macroscopic properties predicted by the models are in good 
agreement with experimental quiet- time observations at 1 A. U. The proton 
velocity distribution function obtained at this radius resembles closely 
that inferred from in- situ proton measurements^ The models may be used 
with increased confidence to predict flow conditions at other heliocentric 
radii presently under experimental investigation or to be explored in 


the near future. 
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I* INTRODUCTION 


I.l The Solar Wind 

The existence of a continuous high-speed outflow of corspuscular 

J 

radiation from the sun, known today as the "solar wind”, is a well 
established fact, first suggested by Biermann ( 1951 ) to explain the 
observed accfeleration of comet tails pointing away from the sun, and 
later predicted by Parker ( 1958 ) as a continuous supersonic expansion 
of the solar corona. 

In his pioneering paper Parker demonstrated that the corona cannot 
exist in a state of hydrostatic equilibrium. Its large extent and elevated 
-temperature, (of the order of a few million degrees), create a pressure 
force distribution that cannot be balanced by the containing effects of 
the sun’s gravitational attraction and interstellar pressure and thus 
expands supersonically into space. This expansion process is analogous 
to the flow of gas through a deLaval nozzle, as pointed out by Clauser 
( 1960 ). 

The existence of the solar wind was verified by the first Mariner 2 
^ result^ (Neugebauer and Snyder, 1962 ; Snyder et al., 1963 ), resolving a 
long standing controversy between the evaporative processes proposed by 
Chamberlain ( 1960 ) or ”solar breeze”, and the hydrodynamic supersonic 
expansion of Parker. The history of the ideas and experimental observations 
that led to the solar wind concept as presently known, has been reviewed 
by Dessler ( 1967 ) and Spreiter and Rizzi ( 1972 ). 

Today, after a decade of rpacecvaft observations, the large scale 
features of the solar wir^d such i \ itf average flow speed, density, 
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composition, electron and proton temperatures and thermal anisotropy 
ratios, are relatively well known. On a smaller scale, it has been 
observed that the coronal expansion is a dynamic process giving rise to 
a multitude of magnetohydrodynamic phenomena such as shock waves, density 
and magnetic field discontinuities, and high-speed stream Interactions, 

In spite of /this dynamic character, a “quiet-state" of the solar wind 
has been associated with low-speed conditions observed at certain times 
to prevail for periods long compared to the expansion time (Hundhausen, 
1972), The observed properties of this "quiet-state" solar wind have 

« 

been summarized in Table I for future reference. 

Reviews of the observational knowledge of the solar wind have been 
given by Ness (1967), Ax.ford (1968) and Hundhausen (1968, 1970), while 
recent measurements concerning transport phenomena, pressure anisotropies 
and other related features, have been reported by Montgomery (1971) and 
Ogilvie et al. (1968). 

The expanding coronal gas is an electrically neutral, highly con- 
ductive plasma and as such it is expected to carry with it the relatively 
weak solar magnetic field. This frozen- in flux combined with solar 
rotation results in the Archimedean spiral structure of the interplanetary 
magnetic field first suggested by Parker and later confirmed by in- situ 
observations by spacecraft in the Venus-Earth-Mars space. Ness and 
Wil ox (1967) reported a particular large scale feature of the inter- 
planetary field; this is its sector structure associated with polarity 
reversals observed during the course of a solar rotation and persisting 
over periods of several solar rotations. 
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TABLE I 

AVERAGE PROPERTIES OF THE LOW-SPEED 
(QUIET-STATE) SOLAR WIND AT 1 A.U. 


Radial Component of Flow Velocity 

Nonradial Component of Flow Velocity 

Proton (electron) Density 

Electron Temperature 

Proton Temperature 

Magnetic Field Intensity 

Solar Ecliptic Longitude of Field 

Proton Thermal Anisotropy 

Total Energy Flux Density 

Electron Heat Conduction Flux Density 


300-325 Km/sec. 

8 Km/sec. 

8,7 cm ^ 

1.5 xlO^OR 
4 xlO^°K 
5 gsmma 
140° 

2 

-2 

0,25 ergs cm sec 

-3 -2 

7 X 10 ergs cm 


1 

-1 

sec 



The spiral configuration of the field results in the transport of 
angular momentum away from the sun thus exerting a retarding torque on 
its outer layers. In addition, a smaller amount of angular momentum is 
transported by the solar wind in the form of an azimuthal velocity 
component in interplanetary space. 

At large distances .from the sun, the momentum flux and magnetic 
pressure associated with the solar wind become comparable to the total 
interstellar pressure. In this region it is expected that the solar wind 
will undergo a supersonic to subsonic transition generating a shock wave, 
^(Axford et al. 1963; Dessler, 1967). 

1*2 Theoretical Models of the Solar Wind 

Since the early work of Parker, (1958, 1960) numerous fluid and 
exospheric models of the coronal expansion have been proposed. Hundhausen 
<1968, 1970,' 1972) has reviewed the general characteristics and conditions 
of applicability for these models, and the accuracy with which they 
predict observed flow conditions at the earth's orbit. 

The as sumption, made in most models that the solar wind behaves 
collectively as a ionized fluid cannot be substantiated in terms of 
if classical plasma theory. The exospheric models of Chamberlain predicted 
very small expansion velocities but later refinements on these models by 
Brandt and Casinelli (1966), Jockers (1970), and Hollweg (1970), produced 
expansion speeds comparable to those obtained from fluid models. 

Nevertheless, other values are in considerable disagreement with observations 
in particular the proton thermal anisotropy ratio and expected behavior 
of He igns in the solar wind. 
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All of the observational evidence indicates a fluid-like behavior 
in the coronal expansion and therefore fluid models are expected to 
give results in better general agreement with observations than exospheric 

I 

(or evaporative) models. 

The general nature of the results obtained from hydrodynamic models 
is the same, that is, supersonic expansion of the coronal gas in inter- 
planetary space. The significant differences among the models result 
from the particular treatment of the energy equation and the inclusion of 
the spiral magnetic field in the analysis. Of particular importance to 
‘the subject of this dissertation are the models of Weber and Davis (1967), 
Urch (1969), Whang (1971a), Wolff et al. (1971) and Whang (1972). 

Weber and Davis developed a one-fluid model with a polytropic 
radial temperature dependence and Included the effects of the frozen- in 
solar magnetic field. The radial expansion velocity is not affected to 
any large extent by the inclusion of the field but a significant retarding 
torque on the sun is predicted as a result of the stress produced by the 

spiral structure of the magnetic field. 

Urch obtained a numerical solution to the magnetohydrodynamlc 
'^one-fluid equations under the assumptions that heat is carried away by 
conduction from a heat source located at the base of the corona and that 
the magnetic field inhibits the transport of thermal energy at tight 
angles to the field. The temperatures predicted by this model at 1 A.U. 
are too high, although other quantities are in agreement with observations. 
The azimuthal velocity at the earth's orbit predicted by these models lies 
in the range of 1-2 K/sec. in disagreement with reported observations 
of 6^10 Km/sec;# 
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Whang considered a radial model of the coronal expansion including 
the spiral magnetic field and showed that under these assumptions, mag- 
netic field expansion energy is continuously converted into kinetic 
energy and thus was able to increase the predicted radial velocity at 
1 A,U* by 17% * Mbdisette (1972) has pointed out that although the magnetic 
energy conversion process described by Whang is indeed operative in the 
solar wind, its effect should not be as large when the azimuthal velocity 
component is taken into consideration in the analysis. We shall consider 
this problem in detail in the first part of this dissertation by in- 
corporating the azimuthal velocity into Whang^s one-fluid model and 
obtain numerical solutions to the resulting system of magnetohydrodynamic 
equations. 

Wolff et al. have proposed that viscosity plays a major role in 
heating up the protons in the solar wind and in this fashion account for 
a non- thermal source required by two- fluid models to obtain reasonable 
proton temperatures at 1 A.U., (Hundhausen, 1970). Although their results 
agree quite well with observations after an empirical function for the 
conductivity coefficient is assumed, the role of viscosity and other 

If 

transport phenomena defined in terms of classical plasma theory is not 
completely understood at the present time. The observational evidence 
points out to the existence of randomizing effects other than Coulomb 
collisions in solar wind but the exact nature of these Interactions is 
not known. 

A theoretical approximation tt> the observed proton velocity distri- 
bution function was obtained by Whang (1971b) and it allows the set of 
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magnetohydrodynamic equations of Chew, Goldberger and Low (1956) to be 
closed relating the proton temperature to the proton heat flux in the 
context of guiding center plasma theory rather than classxcal heat 
conduction. These results were included in a two-fluid, two-region 
anisotropic model (Whang, 1972), capable of providing macroscopic as 
well as microscopic information about the solar wind. 

Weber (1967) and Weber and Davis (1970) have considered the effects 
of thermal anisotropies upon the angular motion of the solar wind under 
certain simplifying assumptions and show that the predicted azimuthal 
velocity at 1 A.U. is five times larger than that predicted by isotropic, 
one- fluid models. The second part of this dissertation will consider the 
solution to the magnetohydrodynamic equations for a two-fluid, t,wO-j.egj.on 
solar wind model which includes Whang’s formulation of the proton velocity 
distribution function to represent proton thermal anisotropy effects 
upon the angular motion. The two-region formulation of the model is 
necessary to avoid the rapid proton cooling problem which has plagued 
most two- fluid models and represents only an approximation to the physical 
processes believed to be responsible for the observed properties of the 


solar wind. 



( 
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II. A ONE-FLUID MAGNETOHYDRODYNAMIC 
MODEL OF THE SOLAR WIND 

II. 1 Basic Assumptions 

In this section we shall obtain numerical solutions to the steady 
state magnetohydro dynamic equations governing the expansion of the 
coronal gas under the assumptions that the solar wind is a perfectly 
conducting fluid, heat is carried away by conduction from the base of 
the corona and that the solar magnetic field depends only on latitude, 
'ignoring its sector structure. 

The assumption that the solar wind behaves as a fluid is based 
upon observed characteristics of the coronal expansion. We shall consider 
this fluid to be inviscid since energy supplied by thermal conduction is 
much larger than that dissipated by viscosity (Parker, 1965) and the 
general form of the viscous stress tensor for low density plasmas in the 
presence of a magnetic field is not well known. Viscous models based 
upon the classical formulation for this tensor have shown greater dis- 
agreement with observations than inviscid models (Whang, Chang and Liu, 
^1966; Scarf and Noble, 1964). The gas will be assumed to be composed of 
fully ionized hydrogen ttith a 5% helium number density content. 

In developing the model we shall closely follow the approach of 
Whang (1971a) , except for the inclusion of the azimuthal component of 
the momentum equation in our analysis. The model will thus represent 
the flow of ionized gas in the sun's equatorial plane and it is further 
assumed that this flow is axially symmetric about the sun’s rotation 


axis. 



MHD Governing Equations for the Model 

The steady state MHD equations of mass, momentum and energy 
conservation, assuming that charge neutrality is maintained in the plasma, 
may be expressed as 

7.Cnu) o (II-2.1) 

(T\nC H ^ (II. 2. 2) 

V. V,W +7q +V. 

+ ^ (ir.2,3) 

where P is Che pressure tensor, N the Poynting vector and X an external 
force function. Other quantities such as the heat flux vector 
magnetic field B, mass density p = mn, temperature T, are represented 
in standard MHD notation and gaussian units will be used throu^out the 
development . 

Maxwell's equations govern the steady state interplanetary magnetic 
field and may be written as 


Vxg 

= o 

(II. 2.4) 

v.n 

ST O 

(II.2.5) 



(II.2.6) 
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V. E s o 

(n.2.7) 


For a perfect conductor in an inertial frame of 

reference the 

electric field is given by 



E - -i- ux B 

(11.2.8) 

and 

the Poynting vector takes the form 


H s. E X 

— 4TT - - 

From (II. 2. 4), (II. 2. 6) and (II.2.8) we obtain 

(II. 2. 9) 


Vx C It >'• §) = 

(II. 2. 10) 

‘and 

1 , JxB = -L C*^xB )>< B 

(II.2.11) 


C -- -- AV^ 


We consider now the spherical coordinate system (r,iy,u;) centered 
at the sun and aligned with the ecliptic plane, shown In Figure la* In 
this system we express the magnetic field B and velocity vector u as 

U = grUf Cr,v) (II.2.12) 

B = firBr CH'P) (h.2.13) 

The mass conservation equation (11.2*1) becomes then 

or 

riUrf^ t: con^r'C. (ii.2.i4) 

since the model is assumed isotropic and one-fluid, the pressure tensor 
is ?iven by 

P =. (2nlVl")V (11,2.15) 

where n is the particle density, k Boltui." s ccnstan’: rr.l I the unit 



Figure 1 


a) The spherical coordinate system. 

b) Definition of the magnetic field angle 0. 
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tensor. It follows that 

7, P vC^nkTj 


introducing (II.2.12) and (II.2.13) in (11.2,10) we obtain 
which implies 

ruh-Buj - ru.M,Br s= 

An analogous procedure with (II. 2. 5) yields 

iz(r^Br) = c, 


(11.2.16) 


(II.2.17) 


(11,2.18) 


(II.2.19) 


or 

r* dr - c**'^*'^^^* (II.2.20) 

The radial and azimuthal components of the magnetic force may be obtained 
from (II. 2. 11) and are given by 

^Q^e)^= ^ (II.2.21) 

and 

t=r (n.2.22) 

r 9r ^ 

where sin ^ has been taken as unity for the equatorial plane model under 


cons iderat ion • 

The external force function X in the case of the solar wind is 


given by the sun’s gravitational attraction 


X !=-er 


-pa 


(II.2.23) 


where G is the universal gravitational constant and the mass of the 
sun. The radial and azimuthal components of the momentum equation thus 

take the forms 
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rUrt^M® 


and 




Ajrr 'dr 


(11,2.24) 


(II. 2.25) 
3 


Here m is the mean mass per particle and n the particle density per cm 
In the case of the energy equation, proceeding in an analogous fashion, 
we obtain 




(II. 2. 27) 


In order to obtain an adequate expression for the heat flux term 
7.^ in (II. 2. 3), we irrust take into account the inhibiting effect of the 
magnetic field upon the transport of thermal energy perpendicular to the 
field lines. Following the approach of Urch (1969), Wolff et al. (1971), 
Whang (1971a), Gentry and Hundhausen (1969), we express the radial 
f component of the heat flux terra as 

where 0 is the angle between the radial direction and the magnetic field 
as shown in Figure lb, and K is the thermal conductivity coefficient. 

Th( energy conservation equation (11.2.3) is thus given by 


I ^ 
•i-irr* 'dr 
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mn 




Ur = o 


(II. 2. 29) 


(II.2.30) 


This equation may be integrated once with the result 

nv^rHSk.- + ^ r*u., 

^ 4TTn r ^ ^ 

=. F 

where F is the total energy flux per steradlan. The second term on the 
left-hand side of (II. 2. 30), not included in Whang's analysis, represents 
the energy flux associated with the rotational motion of the gas. 

The azimuthal component of the momentum equation (II. 2. 25) may be 
integrated directly (Weber and Davis, 1967) to give 

Equations (II.2.18) and (II. 2. 31) may now be used to calculate 

. ^ 4Tmna?c" 4-BfC* (II. 2. 32) 

F 4tr«irtl/r* - 

0uj - S. Uf Br 4 <^ ,) (n.2.33) 

*“ r 4TrmnUr^— 0|.* 

The remaining terms in (II.2.24) are obtained from the above relations 


(II. 2. 31) 


and the equation takes the form ^ « -i 

' — ?5~ 5> ^ ** *>jim«1lrCI-A»> I 


' dlk = ^ 

d.r ^ r 


where 


Uf\ ^ ^ Bu.^ 

4TT»nn0-/l*) 

A'-,:— = ©ry4¥«»nV 


(II. 2. 34) 


(II. 2.35) 

and A is defined as the reciprocal of the radial Alfven Mach number, M^. 

It is Interesting to compare (11.2.34) with the corresponding 
equation obtained by Whang; this will be carried out in Section II. 3. 
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It is convenient to introduce 

+<xn^ « Bu,/Br (n.2.36) 


in equation (11.2,34). Hence 

I ^ r m dr 




dur _ 2 ir 
*■ 1. 


gkT 

m 


(11.2.37) 


The governing equations (II. 2. 30) and (II. 2. 37) may be cast in 


dimensionless form by considering the flow conditions at a particular 
radius. Let us denote the conditions at r = by the subscript indicated. 
,The reason for choosing 3 as the subscript will become apparent as we pro- 
ceed with the development. 

We introduce the following dimensionless variables 

Vs- -y Ws ; e S T/r^ } Zer/rj (11.2.38) 

and dimensionless parameters 7, 5 and 6, defined by 

y*! <SMo/rj^u| y n =. m'Ut|/2kT3 ; (11.2.39) 

Equation (II. 2 . 37 ) may now be written in dimensionless form as follows 


de 


2© 


zigAWW -i 


d? 






(11.2.40) 


where the parameter 

yu, c 43 X Bi^/^irnin3it,r^ (11.2.41) 

is defined as the reciprocal of the radial Alfven Mach mimber at the 
reference radius. 
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The denominator of (II. 2. 40) 

(11.2.42) 

«iU vanish for three sets of values of the independent and dependent 

variables, denoted respectively as 7 (Z 2 ,Vj,e 3 ). 

These three sets correspond to the critical points of (II.2.40), first 
studied by Weber and Davis (1967); each critical point occurs when the 
fluid velocity equals the characteristic propagation speed of a possible 

vave mode iti the tnedium. 

The first critical point, closest to the sun (r=r^). corresponds 
approximately to Parker's critical point where the fluid velocity equals 
the local characteristic thermal speed of the plasma, in our case modified 
by the presence of the magnetic field. The second critical point 
represents the singularity introduced in (II.2.40) by the azimuthal com- 
ponent of the momentum equation. At this point the radial component of 
the fluid velocity equals the local Alfven speed as determined by (II.2.35) 
when = 1. At the third and farthest away from the sun critical point, 
the radial velocity is approximately equal to the local Alfven speed as 
determined by the total magnitude of the magnetic field. Since tan 0 is 
' small in this region, we expect the third critical point to be located 

in the immediate vicinity of the second. 

The general topology of the solution differs little from the one 

given by Weber and Davis and a solution curve extending from the sun's 
surface to large heliocentric distances must pass through all three 
critical points. A schematic representation of this topology is given in 
Figure 2; (V^.Zg) is a node point while (V^,Z^) and (V 3 .Z 3 ) are saddle 
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Figure 2 

Schematic representation of the topology of the solution for 

u near the critical points, 
r 
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points. It is then possible to determine the slope of the solution curve 
at the saddle points from the values assigned to the dimensionless parameters 
Because of this topology, we choose the third critical point as 
the reference radius for the dimensionless equations, Tliis choice will 
eventually determine the success of the numerical integration scheme 
utilized to solve the system of differential equations. Hence, at 2 = 1, 
(11,2.42) takes the form 

I _ s=.o (II. 2. 43) 

from which we obtain for 5 

g ( I-/*) /(^ (11.2,44) 

The denominator of equations (II. 2. 32) and (II. 2. 33) vanishes at the 
second critical point. Since u^ and must remain finite and continuous, 
we require t;hat the numerators must also vanish at this point. Hence 
we must have 

4 TTmn^UrjC^* -i- = o (11.2.45) 

and C* and related by 

c" = (II.2.46) 


vhere the subscript refers to the flow conditions at the Alfvenic crit- 
i 

ical point. 

In a frame of reference rotating with the sun 3 is parallel to u. 


In this frame (Pneuman, 1966) 

Br ' B 


(II.2.47) 


CjJ 


where n denotes the angular velocity of the outer layers, of the sun. 
Introducing (II. 2. 47) in (II. 2. 18) we obtain 


r , 
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(II. 2 . 48 ) 


and 


c*‘= J2r^ 


(II. 2. 49) 


Introducing the dimensionless variables and defining two additional 
parameters a. and where 

<r fz/r^ ; 5* * 

we can write the azimuthal component of the momentum equation as 


(II.2.50) 


iw =. ?Z 4 V4An^ 

where tan 0 is given by 

= ?C<^i2*)/VZC<7*/VZ*) 

At the reference radius (Z = 1), these equations reduce to 


and 


^QrlO/O-/) 


(II.2.51) 

(11. 2. 52) 

(11.2.53) 

(11.2.54) 


We may proceed in analogous fashion with the energy equation (II. 2. 30). 
Introducing the dimensionless parameters 

o( - 2 njVr-r^k/TC cos®54 (11.2.55) 


^ ^ -0 V -/ 

and 


^ di.2.56) 

with the thermal conductivity for ionized hydrogen given by (Spitzer, 


1962) 


we obtain 








(II. 2. 57) 


y 
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(II. 2. 58) 

The dimensionless parameters, 0£, P, 5 h, measure the ratios of 


various energy flows at the reference radius. In addition 

.Va 


" K'^'f/z |5/a /[zkril ^ 


(11*2*59) 


that is, the total energy flux per steradian is proportional to the 
constant K in the thermal conductivity coefficient (Whang, 1971a). In 
reality K is a slowly varying function of the density and temperature 
bf the gas (Braginskii, 1965) but in the case of the solar wind it may be 
assumed constant. 

As shown by Whang, a is not an independent parameter. At the 

dv 

reference radius the numerator of (II. 2. 40) must vanish in order for — 
to remain finite and V continuous across the critical point. Hence, 
from equation (II. 2. 58) we must have 


^ = e< - M - (II.2.60) 

and from equation (II. 2. 40) 

^ (II.2.61) 


The parameter a is thus given by 

^ _ e-yg-tgaVeSg/t-iaaiifs/ci-A) ^ 

and the solutions for V and 0 have the general form 

v= ; O -- - -- 

^ ' ' (II.2.63) 

For convenience and future reference we Siir.nnarlze below the principal 
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equations obtained in this section and the relations among the different 
dimensionless parameters defined. 

Radial Equation of Motion 

g V^- ® - §^V4AnV/2*C»7*^A^*-) 

zimuthal Equation of Motion 

r . 6 a 


dV^ V 
a? " Z 


J 


(II. 2. 64) 


10X1 or nuuj-uLi. 

£ W = ^-2 + ^ +a»i^ 


Ith 


(11.2.65) 

(11. 2. 66) 


lergy Equation 

^a, fi,„2 i\, aVi^->j^S0 ! 

ew v»L (II. 2. ( 

ilations Among Parameters 

I = Cf 70 /C' 7 *A‘‘*V 3 ^ 

_ 2^y^4 ^ g 

S tz 5^9+^/^ H- y 2 

4air^ 5^5 = 5 ^ 


^ J 

(11. 2. 67) 

(11.2.68) 

(11. 2. 69) 


(II. 2. 70) 


cr^ = 14 


(II.2.71) 


Xl^3 The Case of a Purely Radial Expansion 

In Section IIo2 we have obtained for the radial component of the 


momentuni equation 



24 


■3?' “ r 


'4kT <SMo gkr^T , .j -^ , ^^*** 


ttr^- g»^T B3 


(II. 3.1) 


-47r«inC»-A^) 

Whang (1971a), in considering the conversion of magnetic field energy 
into kinetic energy in the solar wind, obtained the corresponding 


equation for the case of a purely radial expansion as 

r^kT <SM© ^ 2kr<dn^ 

r rm 3r 


dttr_ jV 
ar " r 


no 


Ur- 


« -2feT 




u> 


(11.3.2) 


4TTmn - 

We immediately observe that the assumption in (II. 3.1) that 

not reduce this equation to (II. 3.2). The third term in the denominator 

P 

remains divided by the factor (1-A ) which leads to the existence of three 
critical points rather than one, as discussed in II. 2. 

The mathematical source of this discrepancy lies in the assumption 

by Whang that the velocity vector u has the form 

^ (11.3.3) 


while the magnetic field vector is represented by 

B s (II. 3.4) 

In other words, the limiting condition imposed at the onset of 

the analytical development rather than on the final differential equation. 
The limits obtained in each case are different,^ leading to the observed 
discrepancy* 

Physically, the exclusion of u^ from the analysis in the manner 
described above, implies that transverse flo^ perturbations that should 
propagate parallel to the magnetic field as Alfven waves, are ignored. 
Near the sun u »u , A»1 and the third term in the denominator 
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of (II. 3.1) is of opposite sign and smaller magnitude than the corresponding 
term in (II. 3. 2). At large heliocentric distances, A«1 and this term 
reduces to that given by Whang. Consequently, the effects of the magnetic 
field on the flow velocity are not expected to be as large as indicated 
by Whang. Modisette. (1972) has carried out a limited analysis of this 
problem and reached similar conclusions. The numerical results obtained 
in Section II. 4 will show that the effects of the magnetic field on the 
flow are indeed smaller than those predicted by Whang's model. In par- 
ticular the "hose angle" of the interplanetary magnetic field is well 
behaved in the vicinity of the sun, tending to 180 as 


XI. 4 Numerical Solutions for the One-fluid Model 

At the reference radius (Z = 1) equation (II. 2. 64) is of the form 

(— ) and may be evaluated by making use of I'Hopital's rule. The result 
'o 


r T 

s j 

The two solutions of (II. 4.1) correspond to the two possible branches 
of the solution curve for the radial velocity, denoted as and 

'"'in Weber and Davis model and indicated schematically in Figure 2. The 
constants e^^, e^, and f^ are related to the various dimensionless 
parameters previously defined^ as follows 

(I-/*) L 




V 
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(11. 4. 3) 


(II.4.4) 

4 =-%^K'-/^)-']'fil 

(II.4.5) 

where 

d,=.<gc, + 

« 

(II. 4. 6) 

<2^2 ^ q»j"><gC3>f4ny{,4.2SI^^ 

■M&i 

- 

(11.4.7) 

C| s 

(II. 4. 8) 

Ca = ^ 

(II. 4.9) 

f Cg =: 

(II. 4. 10) 

C 4 - 

(II. 4.11) 

^5 = ^^,u 4a _ '2^4a.ix^g^ + 

(II. 4. 12) 
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(II. 4. 13) 




(II. 4. 14) 


(II. 4. 15) 


(II. 4. 16) 


= - •/ Oy*) 

The correct solution of (11.4.1) is given by the (+) sign and corresponds 

to a positive slope at the reference radius. 

It is convenient to express the general solutions for V and 0 as 
'functions of coimnonly used plasma parameters, rather than those given in 
^^■j^ 2.63). For this purpose we Introduce P, the ratio of the thermal 
pressure 2nkT. to the total magnetic field pressure B^/8 tt. Thus at the 

reference radius ^ , 

a *<• 'a « 2 c toS (II. 4. 17) 

^ B’/3T " SA 

arid in general ■ 91^1 

e, ^ © z / VitosVj (11.4.18) 

Making usa of (II.2.68) It follovs that the paranatat 5 will be given by 

•e /'I . ' \+VC I t-L.fi (n.4.19) 

S =Kz-^fgJ-V(T^P!J J 

‘*once and 0^ are specified. The two solutions of (II. 4. 19) represent 
two possible choices for the reference radius. The (+) sign corresponds 
to the outer critical point r=r3 where 0>1. while the (-) sign corresponds 
to the solution that should have been used if we had chosen r=r^, (P=^3^» 
as the reference radius where ^1. Thus the general solutions for 
V and 9 depend on five independent, dimensionless parameters, namely P,, 

Q, H, ^3 and 7. Hence we may write 
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Vs ) i W) 

' * (II. 4. 19) 

From the values assigned to these parameters we can compute the physical 
location of the outer critical point or reference radius, and the flow 
velocity at this point, the results are given by 


^ /SKo?Y’ 

Hr, = ( "SMaQ \<t 

‘ \ rr ^ 


(II.4.20) 


We observe that these quantities are uniquely determined by the particular 
values chosen for the parameters in a given model. Once a solution has 
been obtained, it is not possible to adjust the location of the critical 
point to obtain a best compromise between the predicted temperature and 
flow velocity at 1 A.U. , as the case of strictly radial flow (Whang, 

1971a, 1972). In this sense, we expect the solutions to our system of 
equations to be unique for the particular set of parameters chosen. 

Since the reference radius is located at a certain distance from 
the sun*s surface, two numerical integrations of the differential equations 
starting at the reference radius are required to obtain a complete 
solution: a) An inward (towards the sxm) integration from the outer 
critical point, and b) An outward integration to large heliocentric 
distances. 

The inward integration must pass through two additional critical 
points and presents the most difficult computational problem. Nevertheless, 
the proximity of the Alfvenic critical point to the reference radius 
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considerably simplifies this problem. An inward numerical integration 
starting at the outer critical point will always pass through the 
Alfvenic critical point since this is a node point (Fig. 2) and the 
integration error accumulated over such a short distance is extremely 
small and does not affect the solution in the vicinity of Z^. Thus, 
the inward integration problem reduces to that of finding a solution 
curve starting at the reference radius and passing through the inner 
critical point. This simplification is the most important reason for 
choosing the location of the outer critical point as the reference 
radius. 

Further insight into the behavior of the solutions as a function 
of the values assigned to the parameters in (II. 4. 20), may be obtained 
by considering certain approximations. The five parameters are expected 
to interact to some degree upon the final solution due to the non-linear 
character of the equations. The purpose of the following approximate 
analysis is to find an optimum strategy for the selection of parameter 
values which will satisfy observed conditions in the solar wind, and at 
the same time generate solutions passing through all three critical 
points and satisfying conditions at infinity. 

At the inner critical point (Z^, 0^), the numerator and denom- 

inator of (II. 2. 64) must vanish simultaneously for u^ to remain continuous 
across the point. Hence we must have 

•= O (II. 4.22) 
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and 


(II. 4,23) 


|2e,- = O 

The energy equation (II. 2. 67), at the same point, has the form 

Thus ve only have two independent equations to conipute three unknowns 
and unless we make some assumptions regarding one of the 
unknowns the location of the critical point cannot be determined without 
actually integrating the differential equations. Equations (II. 4 .23) 
through (II. 4. 25) may be normalized to the flow conditions at the inner 
critical point by introducing the following parameters 

<h.4.25) 

r, = = SW.A 

and from (II. 4.23) through (II. 4. 25) we obtain 
wd 


(II. 4.26) 


(II. 4.27) 


(II. 4.28) 


Hence 


These equations are analogous to those obtained for the outer critical 
point and we may write 

_ S.g.yU, /0-^y ai.i.30) 


\ 
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hence, from (II. 4.25) 

0^^ - (II. 4. 31) 

We now make the assumption that in this region the temperature variation 

as a function of radius is given by 

9 (II. 4. 32) 


Hence 


= ml:”'' 


(II. 4. 33) 


This assumption is approximately true in the region considered, for most 
solar wind models (Hundhausen, 1972); we introduce it here for the sole 
purpose of estimating bounds for (||) at the reference radius. Intro- 


ducing (II. 4. 32) in (II. 4. 31), we obtain 




(II. 4. 34) 


Since (jOCcos% ^/a^cos% ^)>1 and must be less than one to represent the 
inner critical point, we must have 


< - 


o,4‘ 


In addition hence 




(II. 4. 35) 


(II. 4. 36) 


^and neglecting the effects of the azimuthal velocity, from (II. 4. 24) 


we may write 


(M), - ^ (11.4.37) 

■^1 

Introducing (II. 4. 32) in this expression, it follows that 

:2 ^ f (II. 4. 38) 

The term in parenthesis may be estimated from previous solar wind models 
and expected conditions at the critical radius and we find that it is 
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less than one. Thus, if Z^<1, we must have uO-1. Equation (II. 4. 35) 
can then be expanded to include this lower bound, 

^ 4.-0.4 (n.4.39) 

At the Alfvenic critical point, equation (II. 2. 64) reduces to 

- - “Z S (IIe4.40) 

Because of the proximity of this point to the reference radius, we may 


write 


-aVWa 


(II. 4. 41) 


(11.4.42) 


Introducing (II. 2. 70), it follox^s that 

“"2 0 

Thus for a given angle 0^, (* 5^)3 almost exclusively on the value 

assigned to the parameter 

It is now possible to formulate an optimum integration procedure 
based on the above results wliich will generate the desired solutions » 

The numerical integration process is optimized by introducing a new 
independent variable X such that 

g"' (II. 4. 43) 

The governing equations may now be written in terms of this new variable 


as follows 


where 


dv „ _ V 
3* ”x 


(11.4,44) 










(II. 4. 45) 
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and 




In addition ^ 

h(x,v,a) 

dx ■ eV* ^ 


where 


’Kx,v,0) -iCv'-OCv+vO^^X-^) + 

2Tx 


^ J-1 + 


and VoB is given by equation (II. 4. 50) below. 


(11. 4. 46) 

(11. 4. 47) 


(II. 4. 48) 


Given an initial set of values for the parameters consistent with 
expected solar wind conditions at the reference radius and such that 
equation (II.4.3S) is satisfied, equations (II. 4. 44) and (II. 4. 47) are 
integrated inwards for X>l+e by means of a fourth-otder, Runge-Kutta 
algorithm. This integration is carried out several times, each time 
adjusting (^)3 by varying Q, such that the numerator and denominator 
of (II. 4. 44) vanish simultaneously or within a small fraction of an 
integration step at the inner critical point. Once this condition has 
been achieved, the integration is allowed to continue towards the sun’s 
surface. The solution thus obtained is valid in the range r^<^<r 3 but not 
beyond. Whang (1971c) has shown that once H is specified, the behavior 
of the solution for 0 for large values of Z is governed by the value 
assigned to the parameter /. We now integrate equations (11.4.44) and 
(II.4.47) in the outward direction (X<l-e) utilizing the value of C 
determined by the inward Integration. Again, the process is repeated 
several times, each time adjusting 7 such that ZO «. The value 
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of 7 thus obtained, is used to find a new value of Q by the inward 
integration procedure previously described and then the ecjuations are 
once more integrated in the outward direction to find the corresponding 
value of 7* 

This iterative procedure is necessary because of the non-linear 
character of the equations^ it is repeated enough times until the inner 
and outer solutions match across the reference radius with typically 
,01% accuracy* The outward integration is then continued to the limits 
imposed by the available computational accuracy and type of computer 
used to perform the calculation. 

For the problems considered in this work, we have utilized Iverson’s 
(1962) APL/360 language because of its unique characteristics, accuracy 
and conversational nature* A listing of the computer programs developed 
to obtain the numerical solution to the one- fluid MHD equations is given 
in Appendix A* 

A problem that has plagued all models that numerically integrate 
the energy equation, is . the extreme accuracy required to specify parameter 
values. This is due to the form of the equations when magnetic field 
inhibited heat conduction is assumed in the analysis and the requirement 
that the solutions must pass through one or more critical points. Our 
model is no exception, although the formulation of the equations in the 
form given by (11,4*44) and (II. 4. 47) was found to reduce the accuracy 
requirements by several orders of magnitude. To obtain a solution that 
passes through all critical points, Q must be determined with typically 
8-digit accuracy; to extend this solution to approximately 3.5 A.U., 7 
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must be determined with 12-digit accuracy, although this figure is 
strongly dependent upon the value assigned to H. If we try to obtain 
numerical solutions beyond this region, the computation time becomes 
prohibitive and other mathematical methods must be considered to obtain 
the desired solutions for the differential equations. 

For large heliocentric distances, it is possible to find approximate 
analytical solutions for our equations in the form of asymptotic series. 

At large r, the direction of conduction heat flow is dominated by the 
.spiral angle of the magnetic field and the conduction heat flux decreases 
much faster than the thermal energy flux. Hundhausen (1971, 1972) and 
Durney (1971) have shown that in this case the flow at large r corresponds 

to an adiabatic expansion with 

0~ (II.4.49) 

On the other hand, the velocity is expected to approach the limiting 
value V=V«; hence in the limit, W, 0-»o and V-*V„. Equation (II. 2. 67) 



(II. 4. 52) , 
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The formal asymptotic expansions of Whang (1972) can now be used to obtain 
a solution for these equations valid for large Z. Thus we write 




(II. 4. 53) 


where e=z”^^^. The leading terms in (II. 4. 47) represent the expected 
behavior of V and 0 as 


and 


It follows that ^ 

de__da’r4+ 

dz “ 3 


(II. 4. 54) 


(II. 4. 55) 


The coefficients are obtained by introducing (II. 4. 53) trough (II. 4. 55) 
in (II. 4. 51). and (II. 4. 52) and setting the coefficient of every power of 
e equal to zero. For the one-fluid model under study we have calcualted 
the first few coefficients as follows 

C„ = y/P ; Cu = -SA/Z5'’ ) =<=> 

C,a = -5C.. CSC,o->f)/zP 

C,J, = -SBCnCSSCa+A^ln^P 

4 

and f 

Ca, =o V CsAo-^Va ;C2j= [MC„CgSC»+A>6«lc,),^A 

‘ (II. 4. 57) 

where 

r,.. 


and 


(II. 4. 58) 

The constant A in (II. 4. 53) is determined from the conditions obtained 
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at the point from the numerical integration such that both solutions 
■join smoothly at this point. Hence 

A ^ (II. 4. 59) 

Once the solutions for V and 0 have been obtained, we can compute 

iw= ^6'^ -t- V-kiA^ (H.4.60) 


(II. 4.61) 

Two solutions to the one-fluid MHD equations have been obtained in 

4 

the manner described above. The parameter values used in each solution 
are given in Table II with the corresponding dependent parameter values 
and flow conditions at the outer critical point. 

We observe that the location of the three critical points in each 
solution is not very different from that given by the polytropic model of 
Weber and Davis. The numerical solutions pass smoothly through all crit- 
ical points and approach the conditions V=V^, W=0 and 0=0 for Z^®. 
Although both solutions give reasonable values for u and T at 1 A.U. 
those corresponding to Solution #2 are in better agreement with quiet 
time solar wind observations ; numerical values for this solution are 
given in Table III for 1.08 £ r/r^ £ 2086. 


II. 5 Discussion of Results and Physical Interpretations 

Figures 3 and 4 show the results obtained fibr u^, u^, T, P, and 
0, for values of the parameters corresponding to Solution #1; Figures 5 
and 6 show the corresponding results obtaiped for Solution #2. 
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TABLE II 

Parameters Solution #1 


H 

c 

P3 

7 

Related constants: 

* a, 


|A 


r 


5 


a 




u 


U03 


(f)3 

(f?3 


0.8 

169.8® 

■ 0.19481 

0.2 

0.10924 

0.46509 
0.96529 
10.035 
0.014885 
0.98384 
263.28 Km/ sec. 
3.91 Km/sec. 
25.237 
24.829 
3.74 

5.01 X lO^OR 
326.83 Km/sec. 
-0.588 
0.167 


Solution #2 

0.85 

169.8° 

0.19531 

0.2 

0.097751 

0.31673 
0.96529 
10.035 
0.01538 
0.98388 
272.98 Km/sec. 
4.198 Km/sec. 
26.234 
25.81 
3.96 

5.39 X 10^®K 
349.86 Km/sec. 
-0.523 


0.173 



1.084 

1.494 

2.065 

2.851 

3.944 

5.465 

7.603 

10.70 

14.99 

21.86 

29.15 

39.16 
52 . 4.7 
71.89 
97.19 

131.2 
181.0 
262.6 

375.3 
499.5 

664.8 

CS 4 .C 

1177 

1567 

2086 


22.29 

47.61 

76.05 

107.2 

138.3 

167.4 

195.2 

220.8 

242.8 

263.9 

277.8 

290.2 

300.8 

310.6 

318.7 

325.5 

331.7 

337.4 

341.6 

344.4 

346.2 

347.4 

348.2 

348.7 
349.1 


2 . 68 x 10 
2.28 " 
1.93 
1.65 ■’ 

1.41 " 


1.19 " 
1.01 ", 
8 . 54 x 10 

7.20 
5.93 


It 


5.10 

4.36 

3,72 

3.12 

2.63 


It 

ti 


- 62.78 

■ 37.70 

- 22.75 

- 13.83 

- 8.54 

- 5.30 

- 3.26 

- 1.97 

- 1.21 

-.700 

-.446 
-.288 
-.187 
-.118 . 
- 7 . 70 x 10 


1.74 

2.18 

It 

- 5.11 

U 

1.36 

1.75 

If 

- 3.36 

If 

1,00 

1.30 

If 

- 2.08 


.735 

9 . 44 x 1 0 "^ 

- 1.36 

ff 

.569 

6.47 


- 8 . 32 x 10 ' 

.435 

4.44 

ir 

- 4.29 

11 

.331 

3.04 

If 

- 2.21 

It 

.251 

2.08 

ft 

- 1.14 

11 

.189 

1.42 

ir 

- 5 . 87 x 10 

.143 

9 . 73 x 10 ^ 

-3.01 

11 


-4 


179.5 

179.3 
179.1 

178.7 

178.3 

177.7 

176.8 

175.6 

173.9 

171.3 

168.7 

165.3 
161.0 

155.2 

148.4 

140.8 

132.1 

122.3 

114.1 
108.7 

104-3 

100.9 

98.2 

96.2 
94.7 


.0189 • 

6101 

1975 

.0162 

1704 

1966 

.0165 

558.8 

1798 

.0190 

208.0 

1546 

.0240 

84.20 

1270 

.0323 

36.25 

1009 

.0454 

16.06 

783.8 

,0669 

7.16 

592.8 

.1000 

3.32 

445.1 

.1593 

1.43 

320.1 

.2280 

.768 

248.2 

.3277 

.407 

191.5 

.4631 

.218 

148.8 

.6518 

.112 

115.0 

.8607 

6 . 02 x 10 “^ 91.74 

1.057 

3.23 

'» 75.48 

1.186 

1.66 

" 63.84 

1.155 

7 * 79 x 10 “'^ 55.66 

.9893 

3.76 

« 51.25 

.7366 

2.10 

" 49.21 

.5309 

1.18 

" 47.99 

/ 

.3747 

6 . 66 x 10 " 47.28 

.2610 

3.75 

" 46.86 

.1804 

2-11 

" 46.61 

.1241 

1.19 

" 46-47 
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Figure 3 

The temperature 5 radial and azimuthal velocities obtained for 
Solution #1 of the one-fluid model, as a function of heliocentric 
distance* 



VELOCITY (ki.i/sec) 



TEMPERATURE {’K) 
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Figure 4 

The plasma f3 and magnetic field angle 0 as a function of helio 
centric distance for Wolution #1, one- fluid model* 


4 



ANGLE ( OECfiECS ) 





Figure 5 


The teinperature, radial and azimuthal velocities obtained for 
Solution of the one- fluid model, as a function of heliocentric 

distance. 
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Figure 6 


The plasma p and magnetic field angle 0 as a function of helio 
centric distance for Solution # 2 , one-' fluid model. 
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The radial component of the expansion velocity is continuous 
across the critical points and increases monotonically from a few tens 
of Km/sec. near the sun, to a few hundred Km/sec. at large helxocentric 
distances. Solution #2 gives 334 Km/sec. at 1 A.U., in good agreement 

with observations during quiet times. 

The azimuthal velocity component first increases with increasing 
distance from the sun's surface due to the tendency of the plasma to 
corotate with the sun. It reaches a maximum around lOr^ and then decreases 
monotonically with increasing distance. The predicted azimuthal velocity 
'at the Earth's orbit for Solution #2 is 1.19 Km/sec., which is of the 
same magnitude as the azimuthal velocities predicted by Weber and Davis, 
Urch and Wolff et al. 

This value of u^ disagrees with reported observations of 6-10 Km/sec 
for the azimuthal speed; nevertheless the uncertainty in these measurements 
is of the same order of magnitude and further work in this area is 

necessary to resolve this conflict# 

Weber and Davis have shown that the characteristic deceleration 
time for the sun due to the torque produced by the magnetic field and 
angular momentum loss can be written as 

r = -3o/%-£2.'t^ 

If we take = -1.25xl0^^gm sec"^ we obtain t=7.9x10^ years as the 
dt 

characteristic deceleration time for Solution #2. As expected this value 

agrees with previous results obtained by other authors* 

Figure 7 shows the variation of the radial Alfven Mach number as a 

2 “2 

function of heliocentric distance. Near the sun A is very large 

V 



Figure 7 


The radial Alfven Mach number ~ / 4TTmn as a function of 
heliocentric distance. Near the sun «1 reducing the effect of the 
magnetic field upon the radial component of the expansion velocity. 




51 


reducing the effect of the magnetic field upon the radial expansion 
velocity as discussed in II. 3. The magnetic field angle 0 given in 
Figures 4 and 6 remains well behaved in the vicinity of the sun and tends 
to 180° for r-r^. A plot of the characteristic Alfven speed B//4 ttp for 
Solution #2 is given in Figure 8 , while Figure 9 shows the radi 
component o£ 

The predicted temperatures at 1 A.O. are l.OtalO^-K for Solution #1 
and 1.54xlO*OK for Solution «2, uhile the plasma p values are .89 and 1.19 
respectively, in good agreement with observations. Table IV summarizes 
the flow conditions predicted by this model at 1 A.U. and for reference 
we have included the predictions of previous one-fluid models as reviewed 

by Hundhausen (1972). 

. .1 £ T — R rtS M and V calculated, are 

The values of u^, u^, T, ^ » 

independent of the value assigned to the constant K in (II.2.57). To 

determine the particle number density, heat flux, magnetic field, kinetic 


and total energy flux, we must assign a value to K. From (II. 2. 55) it 


follows that 




and 



(II. 5.2) 

t ' as 

(II. 5. 3) 

. i/a . ^ f 

0 ^ (4*n-mnyA ) % /2 

(II. 5.4) 


The magnetic and kinetic energy flows per steradian are respectively 

and mn 01 ^ I 5 ^ 5 ) 

^TT ^ ' 

We find that it is not possible to assign a unique value to K that will 
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Figure 8 

The characteristic Alfven velocity V^= as a function 

of heliocentric distance for Solution #2, one-fluid model* 



( km /sec ) 
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• Figure 9 

The radial component of the characteristic Alfven velocity, 
as a function of heliocentric distance for Solution #2, one*- fluid 
model • 



Var (km/sec ) 




TABLE IV 


PREDICTED FLOW CONDITIONS AT 1 A.U. - ONE-FLUID MODELS 



Sol. #1 

Sol. y/2 

Whang 

(1971) 

Urch 

(1969) 

Whang 5t 
Chang (1965) 

Noble & 
Scarf (1963) 

Weber & 
(Davis (1967) 

Radial velocity 
(Km/secO 

317.6 

334.4 

- 302 

371 

260 

352 

~400 

Azimuthal velocity 
(Km/sec 0 

1.02 

1.19 


.623 


— 

1.0 

Temperature 
(deg. K) 

1.04x10^ 

1.54x10^ 

1.5x10^ 

4.39x10^ 

1,6x10^ 

2.77x10^ 

2x10^ 

Magnetic field 
angle (deg.) 

126.2 

127.7 

129.5 

- 

“ 


135° 

Plasma Beta 

.890 

1.19 

1.58 

- 

- 

- 

- 
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give reasonable values for these quantities at 1 A,U. and in the vicinity 
of the sun simultaneously. This would Imply that the assumption K=const. 

Is not valid throughout the region considered. In Figure 10 we show the 
density and magnetic field intensity obtained from Solution #2 for two 
extreme values of K==8xl0 ^ and K=6xl0 ^{ergs cm ^sec deg ). Analogous 
results are given in Figure 11 for the thermal energy flux q* 

The best agreement with observations at 1 A*U. is obtained when 
we choose K==l .0xl0"^(ergs cm ^sec ^deg Table V shows the values 

obtained for these quantities at 1 A,U. for different values of K between 
the two extreme values considered above. 

In Figures 12 and 13 we have plotted the kinetic and magnetic 

energy fluxes per steradian as functions of heliocentric distance for 
-8 “1 -1 -3 5 

K=8xl0 ergs cm sec deg * ♦ It is immediately apparent that only a 
small amount of magnetic field energy is converted into kinetic energy, 
in contrast to the results obtained by Whang, The principal factor 
responsible for the 177e> increase in radial flow speed obtained by Whang 
is the introduction of magnetic field inhibited heat conduction in the 
energy equation. Thermal energy piles up behind the obstruction 
represented by the field, raising the temperature and increasing the 
velocity, (Parker, 1971; Hundhausen, 1972). 

The value of the constant K which gives best agreement with obser- 
vations at 1 A.U, is 0.16 of the classical Spitzer^s value for ionized 
hydrogen. Recent observations of radio-star scintillations indicate that 
the solar wind is highly turbulent, (Jokipii, 1972). These fluctuations 
and waves are an essential aspect of the solar wind and affect the transpoi 



58 


Figure 10 

The particle number density and magnetic field intensity as a 

-8 -7 “1 

function of a heliocentric distance, for K=8xl0 and 6x10 ergs cm 
sec”^deg“^*^. Solution #2, one-fluid model. 
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“ Figure 11 

The thermal energy flux Q, as a function of heliocentric distance 
for K=8 x 10‘® and 6xl0'^ergs cm"^sec‘^deg"^‘^. Solution #2, one-fluid 

model. 



HEAT FLUX (ergs/cm^ sec*) 






TABLE V 


TOTAL ENERGY FLUX, PARTICLE DENSITY. MAGNETIC FIELD 
INTENSITY AND HEAT FLUX PREDICTED BY SOLUTION #2 AT 1 A. U. 


: 10^ - _ 
^sec“^deg"^*^) 

F X 10"^^ 

(ergs sec"^sterad”^) 

-1 

(ergs cm “"see ) 

"-3 

(cm 

B 

(gammas) 

/ *^-2 
(ergs cm “^sec 

6.00 

33.0 

1.487 

35.4 

17.8 

4.62 X lO"^ 

5.00 

27.5 

1.239 

29.5 

16.2 

3x84 X 10 ^ 

A. 00 

22.0 

0.991 

23.6 

14.5 

3.08 X 10*^ 

3.00 

16.5 

0.743 

17.7 

12.6 

2.30 X lO’^ 

1.60 

8.25 

0.371 

8.85 

8.9 

1.15 X 10"^ 

l.AO 

7.70 

0.347 

8.26 

8.6 

1.07 X lO"^ 

1.20 

6.60 

0.297 

7.08 

7.9 

9.25 X 10“^ 

1.00 

5.50 

0.247 

5.90 

7.2 

7.71 X lO"^ 

0.80 

4.40 

0.198 

4.72 

6.5 

6.16 X lO’^ 


e 
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Figure 12 

1 2 2 

The kinetic energy flux per steradian, KEF=^ mnUj.r (u ), 

-8 “1 "1 

function of .heliocentric distance for K=8xl0 ergs cm sec deg 
Solution y/2, one- fluid model. 




as a 
3.5 
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Figure 13 

The magnetic field energy flux per steradian, 

2 2 

2 

“8 —1 -1 — ^ 

as a function of heliocentric distance for K=8xl0 ergs cm sec deg ^ 
The amount of magnetic field energy converted into kinetic energy is a 
small fraction of the total energy flux« 
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coefficients in the plasma reducing the effective heat conductivity along 
the magnetic field lines and increasing the energy exchange rate between 
electrons and protons, (Hollweg, 1972; Hollweg and Jokipii, 1972; Perkins, 
1973). Perkins has argued that Spitzer's conductivity is inapplicable 
in the region where the solar wind becomes collisionless and a reduced 
value should be used instead. The results obtained from the present 
model give support to these hypotheses since the choice of K=6xl0 near 
the sun, corresponding approximately to Spitzer’s value, leads to coronal 
(jejisities in better agreement with observations (see Frgure 10) , while 
the reduced value K=lxl0'^ gives agreement at 1 A.U. where the wind is 
essentially collisionless. These results imply the existence of two 
characteristic regions in the expansion process, with a transitional 
region in between. We shall consider these concepts in greater detail 
in the following chapter and develop a two-region, two- fluid model of 
the solar wind. 


If 


V 
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III. THE EFFECT OF THE PROTON THERMAL ANISOTROPY 
ON THE ANGULAR MOTION OF THE SOLAR WIND 

III.l Introduction to the Problem and Basic Assumptions 

The particle velocity distribution function for a uniform, 
collisionless plasma in equilibrium in the presence of a magnetic field, 
possesses cylindrical symmetry around the field direction and is of the 
form 

•fCS-) = f (iii.i.i) 

vhere C and C denote the intrinsic velocity components parallel and 
II ^ 

perpendicular to the magnetic field* The second moments of (III* 1*1) 
give the parallel and perpendicular pressures and are related to f by 


p,l =r Jmcflfdc 


and 


(III. 1.2) 


Ri, = 

In addition, the parallel and perpendicular temperatures are defined by 

n (<T|i * pH f n px 

and the total plasma temperature is 

‘ T = ( zm 

The third moments of f give the conduction heat fluxes 


(111. 1.3) 

(111. 1.4) 


and 


(III. 1.5) 




mCj^ 


C|| 




wh' ch are identically zero if f is Maxwellian in form. In a frame of 


reference with its direction aligned with the magnetic field, the 
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pressure tensor P will be given by 

" p = pa I -• 

where I denotes the unit tensor. Solar wind observations indicate that 

in the vicinity of the Earth's orbit the proton pressure tensor is anisb- 

trooic with P >P and furthermore q 7 ^ (Hundhausen, 1972) , inq?lying 
II 1 II J- 

that f deviates from the Maxwellian form and the plasma is not in a 
state of thermal equilibrium. Figure lA shows the contour map of a 
typical proton velocity distribution function reconstructed from obser- 
vational data (Hundhausen, 1970). 

The pressure tensor P given by (HI. 1.6) may be used in the 
formulation of a more complete solar wind model if the rates of change of 

P and P are known. Chew, Goldberger and Low (1956) obtained expressions 

II 

for the second moments of the Vlasov equation under the assumptions 


described above which may be written as follows 


fl»K- ''B-' 


and 


(III. 1.7) 


these expressions involve he3.t flux temns the geneml set 

31 

of magnetohydro dynamic equations cannot be closed in terms of known 
moments of the velocity distribution function and thus two additional 
equations are required to determine q^^ and q^. 

Whang (1971d) showed that the proton velocity distribution function 

in the solar wind could be approximated by 


(III. 1.8) 
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Figure 14 

« 

Contour map of the proton velocity distribution function at 1 A.U 
as reconstructed from observational data, (Hundhausen, 1970). The Z-axis 
corresponds to the direction of the magnetic field and points away from 
the sun* 





' 7 / 
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where f° is the bi-Maxwellian distribution function 


(III. 1.9) 


with 


A „ = (ekT.,/™)''* i A, = (2kV ,n) 




and 


(III. 1.10) 

These are the dimensionless forms of the intrinsic velocity components; 
the dimensionless heat fluxes 7^^ and y ^ are defined by 

^T,i s^^ii/Aii (nkT„/a;) ; nlcTl (ui.1.12) 

The function h is an even function of and and in this notation it 
takes the form ^ 

h(?,„sx) ='*^« (|5r-o+2^i (VO 


Using this form of the distribution function it is possible to compute 
the third and fourth moments of the Vlasov equation in terms of lower 
moments, closing the system of MHD equations. Whang has obtained the 
following expressions: 

) ^v\f\^ (III. 1.14) 

» r 3 ± ) - ei . ("tx VD - err.. vtO 

Dt'' mnB - ^ 

Equations (III. 1.7) and (III. 1.14) thus govern the variation of the proton 
temperatures and T^, and proton heat fluxes q^^ and q^ in a collision- 
less heat conducting plasma. 

Near the sun, the energy exchange rate between electrons and protons 
is high and the solar wind behaves as a thermally isotropic one-fluid, 
that, is, the electron temperature equals the proton temperature and the 
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anisotropy ratio is unity. As we proceed away from the sun the plasma 
density decreases and the interaction weakens causing the electron and 
proton temperatures to become different and anisotropic due to the presence 
of the magnetic field. Early models of the solar wind (Sturrock and 
Hartle, 1966; Hartle and Sturrock, 1968) which attempted to incorporate 
this effect by means of classical plasma theory based on binary Coulomb 
collisions failed to predict anisotropy and the observed solar wind 
conditions at 1 A.U. Due to the weak interaction with electrons, the 
protons cool off too rapidly leading to an adiabatic expansion at small 
heliocentric distances. As a consequence the predicted proton temperature 
at 1 A.U. is low while the electron temperature is high leading to values 
of the conduction heat flux much higher than observed. 

Since 'then several mechanisms have been proposed to explain proton 
heating beyond the region in which classical collisions play a dominant 
role. They include collisionless heating by dissipation of hydromagnetic 
waves (Barnes, 1968, 1969; Barnes et al,, 1971; Hung and Barnes, 1973), 
viscosity (Wolff et al., 1971), coronal Alfven waves (Belcher, 1971), MHD 
^pulses (Papadopoulos, 1973) and electrostatic ion cyclotron waves (Toichi, 
1971). Perkins (1973), on the basis of radio-star scintillation observations 
of the turbulent solar wind, has proposed that magnetoacoustic plasma 
Instabilities are responsible for the increased energy exchange rate 
between electrons and protons. At this time the exact nature of the 
Interaction is not known but we may consider an inner region in which 
the proton and electron temperatures are equal and isotropic, and an 
outer region in which the protons become collisionless and their temperature 
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anisotropic. In between there is a transition region in which the protons 
are neither collisionless nor isotropic or one-fluid. Figure 15 helps to 
illustrate these concepts; the dashed lines represent expected solar wind 
conditions as deduced from observations at 1 ^ while the solid lines 

represent results obtained from the two- fluid models indicated* 

The two-region concept has evolved from the theoretical work of 
Hollweg (1970, 1971), Burlaga (1971), Leer and Holzer (1972), Chen et al. 
(1972) and other authors* Whang (1972) has incorporated the two-region 
^approach into a two- fluid model of the solar wind* This model, using the 
proton distribution function (III. 1.7), provides macroscopic as well as 
microscopic information about solar wind protons; the results show good 
agreement with experimental observations. 

The proton thermal anisotropy is of particular importance in the 
study of the solar wind angular momentum since the pressure tensor P 
given by (III. 1.6) will give rise to additional azimuthal forces not con- 
sidered in Chapter II of this work, Weber and Davis (1970) have included 
these forces, as well as viscosity, in their analysis of the azimuthal 
^motion, in an effort to explain the discrepancy between observed and 
predicted azimuthal velocities at 1 A.U. Nevertheless, the form of 
(P^^-P^) used in their calculations was a simple interpolation formula 

f^u-pjL. ^ (II1.1.15) 

in ▼ lich the parameter e is varied between e = 1 at 1 A,U. and e - 0 near 
the sun and P^ represent the total thermal pressure (electron and 
proton) components parallel and perpendicular to the magnetic field* 

Solar wind observations (Montgomery, 1971) reveal that the electron 



Figure 15 


The two regions considered in the two-fluid model of the solar 
wind* The dashed lines represent expected solar wind conditions as 
deduced from observations at 1 A*U., while the solid lines represent 
results obtained from previous two-fluid models* 
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themal anisotropy at 1 A.U. is small, with consequently, the 

form (III. 1.15)) is expected to overestimate electron thermal anisotropy 
effects and lead to larger azimuthal velocities at the earth's orbit. 

In the following sections we shall expand the model developed in 
Chapter 11 to Include the effects of the proton thermal anisotropy in the 
fashion described by Whang. The electrons will be considered isotropic 
throu^out and treated in the same way as in Chapter II. In the inner 
region we will use the equations developed in Sections II. 2 and II. 3 to 
obtain a one-fluid solution for the model up to a transition point where 
the protons will be assumed to become collisionless. From this point out, 
we will use the equations developed below to obtain a two-fluid solution 
in which the proton temperature becomes anisotropic beyond the transition 
point and into the outer region. The general assumptions made in II. 1 
regarding the interplanetary magnetic field, steady state and fluid 
behavior of the solar wind, apply without modification to this case. 

111,2 Governing Equations for the Outer Region 

The equations of mass, momenttun and energy conservation given in 
^‘chapter II have to be expanded to incorporate the anisotropic proton 
pressure tensor P given by (III. 1.6), the second moment equations of 
Chew-Goldberger-Low (III. 1.7) and the third moment equations of Whang 
(III, 1.14). The general forms (II. 2.1), (II. 2. 2) and (II. 2. 3) are valid 
and we need only consider those terms involving pressure, temperature 
and heat flux. 

The velocity u and magnetic field B will be expressed as in (11.2.11) 
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and (II. 2. 12) thus Maxwell's equations apply without modification. Equation 
(III. 1.6) will represent the proton pressure tensor and may be written as 


fr * P-^ I + =|S ^ (Iii.a.i) 

Furthermore, the force coniponent due to the proton pressure is given by 
(Appendix B) 


V. Pp » *^!^<:B.v)B4(^&.v(f^^) 


(III. 2.2) 

The total thermal pressure tensor for our model is then 

P * (nkli + pj.)5 + C (in.2.3) 

vhere T denotes the isotropic electron temperature. The radial and 
e 

azimuthal compotients of the pressure force are thus given by 

(v.P)r = 5 (nk <«)+ SJ + 






r ^ ^ 03. r gP 


(III. 2. 4) 


and 


(V.p) = <^Sj«is:nj^£[nkCTi,.Ta34llii^^=^<:»i}i»'V 

<(+4e.»'^) 4-^ ^ (in.2.5) 

The corresponding components of the equation of motion are obtained by 
introducing (III. 2. 4) and (III. 2. 5) in (11*2.2) and making use of Maxwell* s 
equations. The results are respectively ^ 



79 


(dUu> 


ca^d si rtd 

U*3!»§rl+ HiSl 

\dr ' 

r / 

1 m*M.r ' ^ 

4irmn0lr J J 


and 


(III. 2. 6 ) 


k Chi-ny 

m r ^ 


(III. 2. 7) 


We observe in (III. 2. 7) that the anisotropic proton pressure introduces 
additional terms in the azimuthal part of the momentum equation and it 
cannot be integrated directly as in Chapter II, These terms will tend 
to increase, the total angular momentum in the solar wind. 

We must now find expressions for the terms in the energy equation 
that involve the new form of the pressure tensor, thus we calculate 

V* ^ I ( 111 . 2 . 8 ) 

~i Ol,u ^65 siri^ n ^ C”^ I 

i 

For the heat flux term have 

V-S = -j!5|:r*^r 

^^ 0^0 q Is now composed of three terms: the radie .1 component of the 

electron heat flux and the radial components of the parallel and perpen- 
dicular proton heat flttxes. In a collisionless plasma in the presence 
of a magnetic field, the proton heat flux is 

f p = ■{ C?.i -- ^ 

For the electrons we consider the inhibited conduction heat flux 



(III. 2. 12) 


V I'^Jc * ' ar 

where 0 is defined as in Figure lb, and it follows that 

^ sU 

Ke = KT^ (in.2.13) 

Introducing these results in the energy equation, we obtain after one 
integration 


integration ^ 

‘ (III. 2. 14) 

V CTh-TI)- r'co5^(- 

where F is the total energy flux per steradian. We must now obtain 
expressions for T^, and from the second moment equations of 
Chew-Goldberger-Low and third moment equations of Whang. Under the 
assumptions made for the model under study, these equations take the 
respective forms 


and 


or 




1 - - “2^ 

' “ nkTTi B>/ 

(III. 2. 15) 

_ Br 9 / 

(III.2.16) 




(III.2.17) 

ffiriB^N or -BrJ 

(111,2.18) 


Wltti 
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^ (III. 2 . 19 ) 

-L 3“r j, j.'l . sWW/'auu «w'| 

Bar r ~'UrBr rV* ( ®r ' r 

In Chapter II we have obtained from Maxwell’s equations and conditions 
in the sun's vicinity 

r(Mrd^ - 'Ww hr) 5 = canst (III. 2 . 20 ) 

hence 

dan^ =. X ^ - XZ) (ni.2.21) 


These equations may now be expressed in dimensionless form by introducing 
the dimensionless variables and parameters defined in Chapter II, except 
for the dimensionless temperatures which in this case are given by 


0 « =Te/r^ 6 . ="n.As. ; 6 j. = 'n/ 7 ^ 

(III. 2 . 22 ) 

and dimensionless heat fluxes 



(III. 2 . 23 ) 

The parameter a is now given by 



(III. 2 . 24 ) 


The system of dimensionless equations may be expressed in the compact 


form 


di di ' Z ' di 


(III. 2 . 25 ) 

at, 

Ci7 


where the coefficients are given by the following expressions: 


^ (HI. 2 . 26 ) 

(III. 2 . 27 ) 


ai3, ~ Vcci^/z. \ <3|4 = VsiMy/a (in.2.28) 
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ais = = o •, = v/a. (111.2.29) 

(III. 2 . 30) 

These coefficients (1=1,8) correspond to the radial part of the 

equation of motion. The coefficients for the azimuthal part, ^ 2 j ’ 
are 

SI ^ (III. 2. 31) 

♦ _ 

^12. = - S? [?Y^ 5(©M-^ji)co»2ji^CdS^-'^v/3*3 (III. 2. 32) 

= - ^Vis«ys/Sinjj)/2 (III.2.33) 

t 

^<=*26 =^2.7=0 2.34) 

= (^v/ 22 )C 6 m-©a)s'«^<=*S'V’+ 25 fe WvVz (111.2.35) 

The equations of Chew, Goldberger and Low may be combined with the third 
i 

moment equations of Whang to reduce the number of coefficients to be 
calculated. The modified third moment equations coefficients are used 
in (III. 2. 25) and may be written as 

aai = c 6 %^(\.^Qx%i^^/ 4 |y*) cosV} (m.2.36) 

aji = (m.2.37) 


% 
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= 1 - (III. 2. 38) 

<=>»» =( 2 § 2 )jMi«*ii^li 2 iM _ s:v.y]4SHMSi!isM'V<?» (III. 3. 39) 

<l4, ^ ^ 2 *coJ 5< (j+t.4^)+ (:^) i.«V [(^a<»|t/25vO -'3 aii.2.40) 

+Qa,([l-ea»sV/2Sv'')3 (iii.s.a) 

443 -=0 ; 444 = (0||<2osV/!2^Y*3-J > 44 s =44^ ~ 447 =<3 (III. 2-42) 

44g r= Qjl ^ coS^C^^^-3) -f ( €>tC0S^/2§v') 

(III. 2. 43) 


The coefficients corresponding to the Chew-Goldberger-Low equations are 
given by 


dgl ? 

- ( 2 c oS /v) (c« ^ ir ) 

(III. 2. 44) 

4.52, 

« SfnjzJ (dll- 

(III. 2. 45) 


* J ) ^54 “ 0 

(III.2.46) 

a sc 

- 2H*oSjjf 

... - 


V 
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<^57 = O i 

« 5 8 » - ®ii 4 0 V ) 

061 = Ox-*- Oj.) 

sln^i^^oS^ (^©x -»-H^d5^ Ox^ 

<?65 =0 j 04 j 4 =5 1 5 dfes ff C > 

<l6fc =• ^ 

^&7 s. 

a«8 - - (4 ) C* ■*■ ^^ 4 ) -*- 

e remaining coefficients are associated with the energy 

d are given by . 

sr5/z a 

71 - := ^73 - -s= <^76 ^ ‘> ^77 s= 2 tfoS j* 

JS s C*iV& + 4 ^ + 

^ fix) - 


(111. 2. 47) 

(111. 2. 48) 

(111. 2. 49) 

(111. 2. 50) 

(111. 2. 51) 
equation 

(111. 2. 52) 
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^ <?J.) (III.2.53) 

2 ^ ^ , 

The system of equations of (XII.2.25) has to be solved in the outer region 

of the model to determine the unknown variables V, 9 , 0 ^ 0 » Q Q.- 

® II II ^ 

In the inner region 0 =0 -0 and the equations reduce to those obtained 

e II X 

in Chapter II except for the CGL and third moment equations which were 
not considered. We shall find that the transition point from the inner 
to the outer regions must be chosen such that it lies beyond the three 
critical points previously discussed* Hence, we need not concern ourselves 
with singularities in the differential equations since they occur outside 
of the region for which the system (111*2*25) is assumed valid* 


III *3 Numerical Solutions 

We found in Chapter II that the solutions for V(Z) and 0(Z) depend 
upon five dimensionless parameters ^ given set 

of values assigned to these parameters the solutions in the inner region 
are found in the same manner as in the one- fluid model up to the boundary 
point between the inner and outer regions (Z=Z^^r^/v^), with Z^>1* At the * 

r 

transition point all quantities are continuous and we must specify the 

value of Q and Q in order to obtain a solution to the system (III. 2. 25) 

II J- 

in the outer region. 

We shall assume that the proton heat flux (Q^^+Q^) at the transition 
point is that flxox available at this radius from magnetic field inhibited 
heat conduction by protons as given by the one-fluid solution for the inner 
region* Braginskii (1965) gives for the electron and proton heat conduction 

coefficients respectively 
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S.)6 Tt, jrn^ 

(III. 3.1) 



(III. 3.2) 


where it is assumed that O) T »1 and t , t are the collision times given 

P P e p 


by 


re c: 3.Sx^Q^Tj'^/c>‘Iii>)nc 


Xy> s 3. ox 10 


(III. 3. 3) 
(III. 3. A) 


and furthermore \ , the Coulomb logarithm is taken as 24 * if T^-T^, 


it follows that for the one- fluid model 

TCi = 'K.3( l+o.a4o74) 


(III. 3. 5) 


and 


fl<. = (III. 3. 6) 

From equation (III. 1.12) we obtain for the dimensionless heat fluxes 
and Q 

l _ /i.. 

(111. 3 . 7 ) 

(111. 3 . 8 ) 



Qn = 

2 VH*V ^ 

Oj.= MiTfi' 

V ^2 V ^ 


V 22 

hence at the transition point Z=Z^ 


Clrl = 


Introducing (III. 3 . 5 ) into ( 11 , 2 . 67 ) and requiring that the total thermal 
en ^rgy flux be continuous across we obtain 


(<5„4 (5.)^ = . 078 ^? 64//W ( 






(III. 3 . 10) 


and 
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' y j. y _ .0782*^-1 (iii.sai) 

^ Ids A 

Thus, once a solution for the inner region is obtained we need only 

specify the ratio 7 h at the boundary Z-Z, to obtain the solution in 

tl i 4 

the outer region. It is convenient therefore to introduce an additional 
parameter defined by 

\ = (IH.3.12) 

and the solutions for V, 9^, 0^^, 0^, and take the general form 

-f--f (,^.^3) ^3^ l-J (III.3.13) 

The requirement that the solutions extend from the sun's surface to large 

heliocentric distances with physically meaningful values, restricts our 

freedom to assign arbitrary values to the parameters in (III. 3. 13). As 

discussed in Chapter II. for any given values of H, C 

adjusted to obtain a solution passing through the inner critical point 

while 7 is determined from the condition that the thermal energy flux at 

infinity is assumed zero. Thus we may write (III. 3. 13) as 

-f (III.3.U) 

that is, the values obtained for the quantities in a given model depend 
i 

upon four independent parameters rather than the six previously indicated. 


From (III. 3. 11) and (III. 3. 12) it follows that 

y _ . 073a 8 V 4 / ds N 


(III. 3, 15) 


The numerical integration procedure used to obtain solutions to the system 


(III. 2 .25) is the same as the one described in Chapter II, The system of 


seven simultaneous non-linear differential equations is written in terms 
of the independent variable X=Z ^ and the results are given in Appendix C. 
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For large Z, asymptotic series solutions may be found by introducing 
the formal asymptotic expansions of Whang (1972) in the system (III. 2. 25). 

, / c - 

These forms are given by c. — -n y 

0+ 

©€ = O’*' ^ ^ 

(III. 3. 16) 

Qj. 

<5„ = fc5j^') 

Ox =. 


In addition to the above ve introduce 

\fj= A-,i(uf^7j^'') (III.3.17) 

as the asymptotic form for the azimuthal dimensionless velocity. For 
large Z ve may approximate 

(III.3.18) 

^ V 

W’ en these forms and their derivatives are introduced in (III. 2. 25) and 
the coefficients of every power of e are set equal to zero, we obtain 

Ves 

and the non-zero coefficients for j=l,8 are 


(III. 3. 19) 
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C,4 ^C^A4 Ca-5AzC2.^)Mtf^ i <^«5 =* -SA:^Cz4Mt^ 

^ t<> = C'A^O t" A 4 V »/ ^ 33 

.C,7 ( )/?P) ( - A4 CT47 - ) 

6 ^ 

^.a = O/I^X ) 

(III. 3. 20a) 

Cza =.C '/^0 2C,o(*/A4+gSC,o^25g) 

^23 = +i4SSC|o- fi>y^)-^zCto 3 

<^z4 = Oil 2j g$»C„ /sA*^ 

Czs = (• W[zA. Cw + - 3-«4) + 

2C.O (6?SC,3 4- 3P»C*i- 2/)/,Ci6-v/>a)-4 
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= ( I 7A, C,7 + 2 «> C„ ( + S S ) + 

C.,(3A»4.-2o5SC,3-5^4^» 4SA,C„)-C,4C'^V|J-A4) + 


•AAjCib] 

Cj, c(l/2|Ai)|ll6>q4C,ft.'2C„(?^^^ 4.225SC,4-6A4 Cii + 

qA^Cka)- i.o(44|SC|s+ ArCi4 )-4A,04 + 


4c,s C^^ + ayg)] 

Ct^ = (V-/r)(44^_ 

^l43=-<^io 3 s= - Cii 

C«.= C,.*-C,5+ + 


(III. 3.20b) 


C46 •= 


C„‘-_ <2,s 


= -C„.3AsA4/4gSl/«As- 
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= 

C57 = 

Csa ^ 
C^*> 

067 
063 : 

C73 ; 

C74 


Cg4 ss - Oil 


3 

^ As 


a 52 


-’ Oj^ ^Ojo <+ 2 Og^^ 


-Ci4 - ^5*3 


- C|i - C(g- s. 04 a 

» -20,c- 46^l^v:»/i6 


« -20„ 


= -*ZOio (2C<i*,+ 0,o^_ 2C,3 + A3A4/-<I t^Vco ^6 

(111.3.20c) 

=. (6C,e C,,//)- 20/4 - ) 


: — C&i4/z^ — 2O1C 


\ 
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fii^/zh%^A7 


c, ^ (M^c,,/iVcA,y (6^./^./7)-(sc,oC74 /t; + 


C7^ - (^^Cis/^ 


(III.3.20d) 


where as before 

and ■ (III.3.21) 

s = 

The coefficients are obtained by successive iterations such that the 
ntunerical and asymptotic solutions join smoothly together at a given value 

of Z. 

^ We have obtained two solutions to the system of differential equations 

(111.2.25) in the outer region and the corresponding one- fluid solutions 
for the inner region. The parameter values used in each case are given in 
Table VI, while the computer programs developed for the two-region model 
ar^ given in Appendix D. The one- fluid solutions are obtained by the same 
procedure described in Chapter II. The transition point from the one-fluid 
formulation to the two- fluid description is chosen at Z=3.3 (-0.4 A. U.) 
and the ratio of 7^^ to 7^ (H) as 4.62 in both cases. These values were 
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Parameters 

TABLE VI 
Solution #3 

Solution #4 

H 

0.8 

0.85 

^3 

169.8° 

169.8° 


.1948 

.1953 

^3 

.200 

.200 

y 

.10911 

.097257 

n 

4.62 

4.62 

Related constants 

A 

a 

.4645 

.31455 


.96529 

.96529 

\ 

10.035 

10.035 

6 

.01488 

.01538 

a 

.98384 

.98388 

^1- 

263.38 Km/sec. 

273.44 Km/sec, 

“ 0)3 

3.92 Km/sec. 

4.20 Km/sec, 

^3 

25.247 

26.278 Tq 

CVJ 

u 

24.839 

25.854 Tq 

^1 

3.75 

3.965 

^3 

5 . 01 x 10 ^ 

5.41x10^ 

u 

326.95 Km/sec. 

350.45 Km. sec 

CO 



<i>3 

-.58753 

-.5184 

<f)3 

.16748 

.17314 


16748 


17314 
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TABLE VI (CONTINUED) 


Parameters 


Solution #3 


Solution #4 



421 

6,1625 

17.346 

.7702 

3.6174 

.03617 

2.6172 


752 

12.035 

20.067 

.8056 

6.9959 

.07407 

2.5661 
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selected to obtain reasonable agreement with experimental observations 
for the proton temperature and anisotropy ratio at 1 A.U. 

III. 4 Results and Physical Interpretations 

The results obtained for the radial velocity u^, azimuthal velocity 
electron temperature T^, and proton temperatures and are 

shown in Figure 16 for Solution #3 and Figure 17 for Solution #4, The 
radial velocity solutions are essentially the same obtained previously 
for the one-fluid models since similar values of the parameters have been 
used in the calculations. 

The azimuthal velocity solutions for the inner region are of the 
same general form as in the one- fluid models; in the outer region the 
effect of the proton thermal anisotropy is to increase the azimuthal 
speed as shown in the figures. The dashed curves represent the more likely 
physical situation rather than the abrupt transition predicted by the model. 
Since the fluid is assumed invlscid, the increase in azimuthal velocity is 
due solely to proton thermal anisotropy effects. The predicted azimuthal 
velocities at 1 A.U. are 1,44 Km/sec. for Solution #3 and 1.68 Km/sec. 
for Solution #4. These values should be compared with those obtained from 
the one-fluid models in Chapter II, 1,02 Kra/sec. and 1.19 Km/sec. 
respectively. Thus the Increase in azimuthal velocity at 1 A.U. is of 
the order of 0.5 Km/sec. 

Weber and Davis (1970), in considering the effects of thermal 
anisotropy and viscosity in the solar wind, incorporated an ad-hoc 
relationship for P which tends to overemphasize the effects of the electron 
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Figure 16 

The radial and azimuthal velocities, the electron temperature 
and parallel and perpendicular proton teTJiperatures obtained for Solution 
#3 of the two-fluid model, as a function of heliocentric distance. The 
dashed lines represent a possible physical situation for the azimuthal 
velocity in the transitional region. 



VELOCITY ( km/sec) 



TEMPERATURE (*K) 



Figure 17 

The radial and azimuthal velocities, the electron temperature 
and parallel and perpendicular proton temperatures obtained for Solution 
#4 of the two-fluid model, as a function of heliocentric distance. The 
dashed curve represents a possible physical situation for the azimuthal 
velocity in the transitional region. 


4 



VELOCITY (km/sec) 



TEMPERATURE (“K) 



anisotropy upon the angular motion. They obtained Km/sec. at 1 A.U. 

although this was accomplished assuming a cne-fluid model in which the 
viscous stress is greatly enhanced by the elevated proton temperature 
assumed. 

More recently Urch (1972) has given a perturbation solution to the 
one~ fluid, isotropic MHD equations, which predicts a mean azimuthal 
velocity of 1-2 Kra/sec at 1 A.U., with excursions of ~ +10 Km/sec. caused 
by the rotating sector structure of the magnetic field. Since ours is a 
steady state model, we cannot calculate time -dependent effects but the 
range of azimuthal velocities that can be considered at the reference 
radius for which physically meaningful solutions can be obtained, is 
considerably smaller than the excursions indicated by Urch. 

We can conclude from the results obtained in this Chapter, that the 
effects of thermal anisotropies upon the angular motion of the solar wind 
are relatively small and cannot increase the predicted azimuthal speed at 
1 A.U. to values in agreement with observations. 

The temperature profiles obtained for the inner region are analogous 

to those calculated for the one- fluid models. In the outer region where 

*the proton thermal anisotropy is allowed to develop, the ratio of to 

T increases rapidly with increasing heliocentric distance reaching a 
pi 

maximum value of ~1.7 at 200 solar radii for both solutions. This ratio 
then decreases monotonically and becomes less than one for large heliocentric 
distances, as shown in Figures 18 and 19. The total proton temperature is 


given by 


=• -4- lp„^/2> 


(III. 4.1) 


/ 
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Figure 18 

The proton thermal anisotropy ratio and the proton to 

electron temperature ratio T^/T^ as a function of heliocentric distance 
predicted by Solution #3, two-fluid model. 






•X9/ 
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Figure 19 

The proton thermal anisotropy ratio proton 

electron temperature ratio T^/T^ as a function of heliocentric distance 
predicted by Solution #4, two- fluid model. 


V 




fro/ 
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The dashed lines show the probable physical situation rather than the 

idealization assumed in the model. 

It is of interest to consider the possibility that wave inodes 
associated with plasma instabilities may be excited in the plasma due to 
the anisotropic proton pressure. Two particular types of instabilities 
warrant consideration: the firehose and mirror instabilities. These 

instabilities will occur if the following criteria are satisfied, 
(Clemmow and Dougherty, 1969; Krall and Trivelpiece, 1973) 



1 > (firehose) 

(III. 4. 2) 

Ipll 

) > 0/lPr^) (mirror) 

(III. 4. 3) 

vhere 


(III. 4. 4) 


and we have neglected the effect of the electrons since they are assumed 
Isotropic. When the instability criteria given by these equations are 
in^osed on our solution we find that (III. 4. 2) and (III. 4. 3) are nowhere 
satisfied and hence no instabilities are expected to occur in the plasma; 
the magnetic field pressure is everywhere greater than the proton thermal 


4 


pressure. 

The plasma P 
and 21, where 


and magnetic field angle 0 are .given in Figures 20 

p <in.4.5) 


As in the case of T and T , the P curve shows an abrupt slope change 

Pll P-*- 

at the boundary between the inner and outer regions caused by the 


idealizations assumed in the model. 
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Figure 20 

The plasma 3 and magnetic field angle 0 as a function of helio- 
centric distance. Solution y/3, two-fluid model. 



AN6L : 
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Figure 21 

The plasma p and magnetic field angle 0 as a function of helio~ 
centric distance. Solution #4, two-fluid model. 
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Flow conditions predicted by the two-region model at 1 A.U. are 
summarized in Table VII and compared with results obtained from previous 
two- fluid models. The quantities show generally good agreement with 
quiet-time solar wind observations. Tables VIII and IX give detailed 
numerical results for Solutions #3 and #4 in the inner and outer regions. 
These two solutions represent typical limits of observed electron temperature 
conditions at 1 A.U. 

The quantities obtained and given in the tables, are independent of 
the value assigned to the thermal conductivity constant K, (see Chapter II). 
To compute the density, energy flow and magnetic field intensity we have 
chosen K = l.OxlO'^ for Solution #3 and K = 1.0x10 ^ for Solution #4; 
these values of K give results that agree reasonably well with experimental 

observations at 1 A.U. 

Figures 22 through 25 show the magnetic field intensity, density 
and heat fluxes predicted by the present model as a function of helio- 
centric distance. Table X summarizes the values predicted for these 
quantities at 1 A.U. For completeness. Figures 26 through 29 show the 
^magnetic field and kinetic energy flows per steradian obtained in each case. 
As in Chapter II, we find that the amount of magnetic field energy converted 
into kinetic energy is small and has little effect upon the final expansion 
velocity. 

The microscopic properties of the solution for the proton distribution 
function are determined by the values obtained for 7 and 
From equations (III. 3. 7) and (III.3.8) it follows that 
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TABLE VII 


PREDICTED FLOW CONDITIONS AT 1 A.U. - WO FLUID MODELS 


y 

Sol. #3 

Sol #4 

Whang 

(1972) 

Wolff et al., 
(1971) . 

Hartle & 
Sturrock, (1968) 

Radial Velocity 
(Km/sec.) 

317.7 

335.2 

331 

303 

250 

Azimuthal Velocity 

1.44 

1.68 

- 

1.8 

- 

(Km/sec. ) 

Electron Temperature 
(°K) 

1.34x10^ 

1.98x10^ 

1.52x10^ 

2.03x10^ 

3.5x10^ 

Parallel Proton 

9.41x10^ 

1.17x10^ 

1.19x10^ 

- 

- 

Temperature, (®K) 

Perpendicular Proton 

6.41x10^ 

6.89x10^ 

5.30x10^ 

- 

- 

Temperature, (®K) 

Total Proton Temperature 

6.75x10^ 

8.5x10^ 

7.5x10^ 

4.0x10^ 

4.4x10^ 


Proton Thermal Anisotropy 

1.73 

1.70 

2.23 

- 

• 

Ratio 

Magnetic Field Angle 

123.3° 

127.7° 

130.1° 

125° 

- 


Plasma P 


860 


1.09 


1.09 



TABLE VIII 


NUMERICAL SOLUTION FOR THE INNER REGION 
SOLUTION #3 


vlr^ 

Ur 

(Km/ sec. ) 

UO) 

(Km/sec, ) 

T 

(°K) 

dQ 

dZ 

0 

(deg.) 


1.274 

36.91 

2.22 

2.51x10^ 

-52.72 

179.4 

.01879 

1.843 

66.29 

2.80 

2.06 " 

-29.54 

179.2 

.01801 

2.626 

98.64 

3.34 

1.71 

-17.15 

178.8 

.02040 

3.733 

129.5 

3.88 

1.42 " 

-10.21 

178.3 

.02612 

5.459 

165.7 

4.21 

1.17 " 

-5.699 

177.6 

.03573 

8.240 

198.2 

4.47 

9.39x10^ 

-3.092 

176.4 

.05454 

12.00 

223.8 

4.49 

7.65 " 

-1.772 

174.9 

.08318 

19.08 

250.0 

4.22 

5,90 " 

-0.892 

172.1 

.1434 

28.05 

267.9 

3.78 

4.71 " 

-0.500 

168.7 

.2265 

37.68 

279.3 

3.35 

3.94 " 

-0.321 

165.3 

,3192 

50.50 

288.9 

2.89 

3.29 '* 

-0.205 

160.9 

.4412 

69.18 

297.6 

2.40 

2.70 " 

-0.121 

155,0 

.6073 

84.18 

302.2 

2.10 

2.40 •' 

-0.079 

150,7 

. 7301 



% 




t' 
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TABLE VIII (CONTINUED) 

NUMERICAL SOLUTION FOR THE OUTER REGION 
SOLUTION #3 




'r© 

"r 

(Kin/sec) 

(Km/ sec) 

T 

Pll 

(OK) 

>.61 

302.6 

2.13 

2.37x10^ 

L7.4 

309.1 

2.25 

1.87 " 

52.9 

314.3 

1.88 

1.33 •• 

29.7 

318.4 

1.32 

8.46x10 

15.9 

321.2 

0.89 

5.13 " 

63.6 

323.3 

0.538 

2.64 " 

17.0 

324.6 

0.379 . 

1.57 '• 

12.3 

325.4 

0.278 

9.10x10^ 

093 

325.9 

0.208 

5.21 " 

455 

326.2 

0.157 

2.96 " 

936 

326.5 

0.119 

1.68 •' 

577 

326.6 

0.090 

9.53x10 

431 

326.7 

0.068 

5.38 " 

566 

326.8 

0.052 

1.96 '* 

078 

326.8 

0.039 

1.71 " 


T 

pi 

(°K) 

2.33x10^ 
1.34 " 
7.89x10^ 
4.89 ** 

3.25 " 

2.07 " 

1.51 •• 

1.12 

8,33x10' 

6.22 " 




T 

dee 


Q, 

e 

dZ 

It 


(OK) 

1.29x10“^ 

5.41x10"^ 

2.38x10^ 

-7.63x10"' 

5.51xl0’ 

1.31 ” . 

1.99 " 

-4.91 ” 

1 75 ’• 

4. 38x10"^ 

2.89x10"^ 

6.13x10" 

1.62 ” 
1.27 " 

-3.23 '* 

-2.12 •• 

1.06 " 

1.54 ” 

9.61x10^ 

-1.41 ” 

1.66x10 y 

3.15xl0"g 

9.96x10 

6.33 '• 

-9.71x10“ 

3.65x10"' 

8.32x10" 

4.24 •• 

-4.78 ” 

3.18 " 

2.86 ” 

-2.39 " 

1.93 ” g 

1.02 . 

1.94 " 

-1.20 ” 

4.57x10" 

3.29x10 

1.32 ” 

-6.14x10" 


-4 


4.66 

3.49 

2.62 

1.48 

1.47 


•• 1.08 

»• 2.59x10 

" 6.20x10 

■' 1.48 

" 3.55x10' 


-10 

-11 


•12 


1.05 
3.39x10 
1.08 •» 
3.48x10 
1.11 


-10 


-11 


9.00x10' 

6.14 •• 
4.19 ’* 

2.86 •• 

1.95 " 


-3.13 
- 1.60 
-8.20x10 

-4.2 •• 

-2.15 ’* 


-5 


0 

(deg) 


150.4 

142.5 

133.6 
124.3 

116.6 

108.9 

104.5 

101.0 

98.34 

96.29 

94.74 

93.56 

92.68 

92.01 

91.51 


P 


.736 

.840 

.891 

.843 

,715 

.512 

.362 

.252 

.175 

.121 

.084 

.059 

.041 

.029 

.020 



TABLE IX 


NUMERICAL SOLUTION FOR THE INNER REGION 
SOLUTION #4 



u 

r 

(Km/ sec) 

u 

U) 

(Km/sec) 

T 

(°K) 

is 

dZ 

0 

(deg.) 

P 

1.206 

31.88 

2.16 

2.5410^ 

-52.95 

179.5 

0.0176 

1 ,665 

56.61 

2.68 

2.15 

It 

-31.81 

179.3 

0.1600 

2.301 

86.32 

3.20 

1.83 

II 

-19.25 

179.0 

0.0171 

3.179 

117.8 

3.68 

1.56 

II 

-11.70 

178.6 

0.0204 

4.409 

147.3 

4.17 

1.32 

It 

- 7.30 

178.1 

0.0266 

6.100 

177,0 

4.49 

1.13 

II 

C 

- 4.50 

177.4 

0.0360 

8.462 

203.6 

4,71 

9.61x10-^ 

- 2.78 

176,5 

0.0511 

11.92 

228.2 

4.77 

8.08 

II 

- 1.69 

175.1 

0.0759 

16.92 

250.0 

4.63 

6.76 

tl 

- 1.01 

173.2 

0.1161 

24.98 

270.6 

4.31 

5.53 

It 

-0.577 

170.2 

0.1880 

32.24 

282.6 

3.92 

4.86 

It 

-0.380 

167.7 

0.2580 

43.44 

294.6 

3.46 

4.16 

II 

-0.240 

163.9 

0 . 3720 

58.40 

304.8 

2.97 

3,57 

te 

-0.149 

159,2 

0.5280 

79.65 

313.9 

2.46 

3.07 

II 

-0.085 

153.0 

0.7470 

87.62 

316.4 

2.32 

2.95 

ft 

-0.069 

150.9 

0.8278 
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TABLE IX (CONTINUED) 


NUMERICAL SOLUTION FOR THE OUTER REGION 
SOLUTION #4 



u 

(Km/sec) 

u 

(Km/Wec) 

T 

PH 

(OK) 

T 

PJ. 

(°K) 



(°K) 

dee 

dZ 

0 

(deg) 

P 

89.10 

316.9 

2.35 

2.92x10^ 

2.86x10^ 

2.28x10"^ 

9.57x10""^ 

2.93x10^ ■ 

-6.61x10"^ 

150.5 

.835 

122.2 

325.0 

2.52 

2.26 " 

1.63 " 

9.22x10"^ 

2.38 ” 

2.56 ” 

-4.28 " 

142.8 

.975 

169.6 

331.5 

2.10 

1.58 •• 

9.66x10^ 

2.92 " . 

5.50x10" 

2.22 " 

-2.85 ” 

134.0 

1.07 

239.1 

336.8 

1.48 

9.96x10 

5.88 " 

7.40x10" 

1 23 ” 

3!20x10"® 

1.87 " 

-1.95 ” 

124.8 

1.08 

328.8 

340.5 

1.01 

6.00 •' 

3.89 •• 

1.83 '• 

1.55 " 

-1.41 " 

117.0 

.997 

526.6 

344.4 

.559 

2.60 " 

2.24 •• 

2.03x10"^ 

4.64x10"^ 

1.07 ”, 

-8.42x10"^ 

107.8 

.749 

752.9 

346.6 

.376 

1.32 •• 

1.52 •• 

3.62x10"^ 

8.34x10" 

1.09 ” . 

3.49x10" 

7.42x10 

-5.89 " 

102.7 

.535 

1002 

347.8 

.280 

7.60x10 

1.12 ”, 

5.01 " 

-3.31 " 

99.69 

.372 

1333 

1775 

348.7 

349.2 

.211 

.160 

4.33 ” 

2.46 

8.38x10^ 
6.26 ” 

1.95 •• 

4.62x10 

1.11 ”.g 

3.55x10 

3.40 " 

2.31 " 

-1.67 " , 

-8.51x10" 

97.33 

95.52 

.257 

.178 

2363 

3145 

4186 

349.6 

349.9 

350.0 

.121 

.092 

.070 

1.39 » 

7.88x10 
4.45 •' 

4.69 ” 

3.51 " 

2.64 ” 

1.10 " - 
2. 62x10" 
6.27x10" 

1 13 ** 

-10 

3.62x10 
1 15 

• -11 
3.69x10 

1.57 " 

1.07 '• 

7.32x10 

-4.34 ” 
-2.21 " 
-1.13 ” c 

94.16 

93.13 

92.35 

.123 

.085 

.059 

5572 

350.2 

.053 

2.51 " 

1.98 ” 

1.50 ” 
3.59x10 

4.99 " 

-5.81xl0" 

91.77 

.040 

7416 

350.3 

.040 

1.42 " 

1.48 ” 

1.17 " 

3.41 " 

-2.98 " 

91.33 

. .028 
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Figure 22 

The particle number density and magnetic field intensity as a 

, „ ,^-7 -1 -1. -3.: 

function of heliocentric distance, for K=1.0xl0 ergs cm sec deg 
Solution #3,' two-fluid model. 





B (GAMMAS 
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Figure 23 

The electron conduction heat flux and proton heat fluxes 

-7 -1 

and Q as a function of heliocentric distance, for K=1.0xl0 ergs cm 
j. 

sec"^deg"^'^. Solution #3, two-fluid model. 





Q (ergs cm sec 
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' Figure 24 

The particle number density and magnetic field intensity as a 

-7 -1 -Ij -3.5 

function of heliocentric distance, for K=1.0xl0 ergs cm sec deg 
% 

Solution #4, two-fluid model. 


h' 




N ( cm*^) 
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Figure 25 


The electron conduction heat flux end proton heat fluxes 
0 as a function of heliocentric distance, for K-l.OxlO ergs cm 
Solution #4, two-fluid model. 




TABLE X 

TOTAL ENERGY FLUX, PARTICLE NUMBER DENSITY, MAGNETIC FIELD INTENSITY AND 

HEAT FLUXES PREDICTED BY THE SOLUTIONS AT 1 A.U, FOR K = 1 x lo"^ (ergs- 

. -1 -1 , -3.5. 

cm -sec -deg ). 



Solution #3 

Solution #4 

F X 

(ergs-sec -sterad ) 

5.68 

5.54 

F / 

^ -2 -1 
(ergs-cm -sec ) 

0.256 

0.249 

(cm 

7.31 

5.85 

B 

(gammas) 

7.69 

7.24 

-2 -1 

(ergs-cm -sec ) 

4.26x10"^ 

l.lOxlO'^ 

/ II -2 - 1 , 

(ergs-cm -sec ) 

7. 00x10*^ 

1.15x10"^ 

q 

X -2 -1 

(ergs-cm -sec ) 

l.OOxlo'^ 

1.94xl0"^ 



Figure 26 


The kinetic energy flux per steradian, KEF - 

-7 

function of heliocentric distance, for K=l*0xl0 ergs 
Solution #3, two- fluid model. 
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Figure 27 


The magnetic field energy flvtx per steradian, 
2 2 

MEF - — ; — ( u sin 0 - u sin 0 cos 0 ) 
4tt r tu 

as a function of heliocentric distance ^ for K=1.0xl0 
Solution #3, two-fluid model. 
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Figure 28 

The kinetic energy flux per steradian as a function of helio- 
centric distance for K=1.0xl0^^ergs cm'^sec ^deg Solution #4, 

two- flu id model. 
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Figure 29 

The magnetic field energy flux per steradian as a function of 
heliocentric distance for K=leOxlO ergs cm sec deg * Solution #4, 
two- fluid model. 
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Ym ^ ) (I1I.4.6) 

^^x•=S^ (ye*7§/e») (III.4.7) 

&± 

The values of 7 , 7 , T and T obtained for various selected heliocentric 
II -L Pll P-L 

distances are given in Table XI and Figures 30 through 43 show plots of 
constant contour maps of the proton distribution function for the parameter 
values given in the table* The velocity scale for each map has been 
normalized to the local characteristic thermal velocity as determined by 
the perpendicular temperature. The axis OZ is parallel to the magnetic 
field direction, facing outwards from the sun* 

The maximum value of the distribution function is attained at point 
0 and the triangle denotes the point in velocity space where the proton 
intrinsic velocity is zero. A comparison of the contour maps obtained 
at 1 A*U, with that given in Figure 14 shows that the form of the distri- 
bution function used in the analysis can adequately represent observed 
solar wind properties. 

^ Whang has given several scale times obtained from the purely radial 

model such as deflection time, equipartition time and expansion time. 

Since the general features of the solution affecting the calculation of 
these times are not very different from those obtained by Whang, we shall 
not repeat the computation here. The equipartition time between electrons 
and protons is much larger than the expansion scale time implying that 
thermal equilibrium between the two fluids cannot be maintained by Coulomb 
interactions alone* We must resort to the physical phenomena mentioned 
previously to explain the increased proton heating within the solar envelope. 





TABLE XI 


PARAMETER VALUES WHICH DETERMINE THE PROTON DISTRIBUTION FUNCTION 
AT SELECTED HELIOCENTRIC DISTANCES 



107 (.5 A.U,) 

.378 

.067 

2.02x10^ 

1.57x10^ 

1.28 

.532 

.102 

2.54x10^ 

2.06x10^ 

A 

1.23 

214 (1 A.U.) 

.403 

.0498 

9.40x10^ 

5.42x10^ 

1.73 

.601 

.0867 

1.17x10^ 

/ 

6.89x10 

A 

1.70 

428 (2 A.U.) 

.384 

.0437 

3.05x10^ 

2.28x10'^ 

1.33 

.616 

.0819 

3.83x10 

2.86x10 

1.33 

1112 (5.2 A.U.), 
Jupiter 

.269 

.0444 

5.21x10^ 

8.33x10^ 

.625 

.578 

.0831 

6.2x10^ 

1.01x10^ 

.613 

2033 (9.5 A.U.), 
Saturn 

.258 

.0455 

1.54x10^ 

4.46x10^ 

.345 

.553 

.0842 

1.89x10^ 

5.48x10^ 

.345 

4066 (19 A.U.), 
Uranus 

,256 

.0463 

3x85x10® 

2 .22x10^ 

.173 

.346 

.085 

4.75x10® 

2.73x10^ 

.174 

6420 (30 A.U.), 
Neptune 

.256 

.0466 

1.55x10^ 

1.4x10^ 

.110 

.546 

.0853 

1.90x10^ 

1.73x10^ 

.110 
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Figures 30 through 36 

The proton velocity distribution function predicted by Solution 
#3 for selected heliocentric distances. The velocity scale has been 
normalized to the local characteristic thermal velocity determined by 
the perpendicular temperature. 
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Figures 37 through 43 


The proton velocity distribution function predicted by Solution 
#4 for selected heliocentric distances. The velocity scale has been 
normalized to the local characteristic thermal velocity determined by 
the perpendicular temperature. 
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IV. SUMMARY AND CONCLUSIONS 

This dissertation has considered the solution of the steady state 

magnet oh ydrodynamic equations governing the supersonic expansion of the 

solar corona into interplanetary space under various assumptions regarding. 

the form in which proton thermal energy is carried away from the sun. 

The fluid has been assumed to be inviscid and the flow axially 

symrnetric about the suxi^s rotation axis* 

In Chapter II we have obtained detailed numerical solutions to the 

* one-fluid formulation of the MHD equations under the assumption that thermal 

energy is carried away by heat conduction from a thin shell heat source 

located at the base of the corona. The effects of the angular motion of 

the solar wind are included in the model as well as a complete description 

of the magnetic field, leading to the existence of three critical points 

through which the solution must pass in order to extend from the sun's 

surface to large heliocentric distances. The magnetic field is further 

assumed to inhibit the flux of thermal energy perpendicular to the field 

-4/3 

lines, leading to an adiabatic expansion at large r with T ~ r 

The values predicted for the flow quantities at 1 A.U. are in good 
agreement with quiet-time solar wind observations except for the azimuthal 
component of the expansion velocity which is approximately a factor of 
five smaller than indicated by reported observations. This discrepancy 
may be due in part to the large degree of uncertainty associated with the 
experimental values and further work in this area seems to be necessary. 

A two- fluid formulation of the MHD equations was obtained in Chapter 
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XII in which the protons are assumed to become collisionless and 
anisotropic beyond an arbitrarily selected radius; the evolution of the 
proton temperature and heat flux can then be described by the Chew- 
Go Idberger-Low theory and third moment equations of Whang, leading to a 
closed set of differential equations which admits numerical and asymptotic 
solutions. These equations were then applied to a solar wind ; model 
consisting of two regions: a) An inner region in which the energy exchange 

rate between protons and electrons is sufficiently high such that their 
temperatures are essentially equal and isotropxc; under these conditions 
the model is adequately described by the one-fluid formulation of the MHD 
equations, and b) An outer region in which the protons are assumed to 
become collisionless and anisotropic beyond a given radius. The electrons 
are assumed .everywhere isotropic and the associated heat flux due to 
conduction alone. The two-fluid formulation of the MHD equations is 
utilized in the outer region to obtain numerical and asymptotic solutions 
for the flow quantities throughout interplanetary space. In addition, 
the formulation of the CGL-Whang moment equations allows us to obtain 
microscopic information about the proton distribution function for various 
heliocentric distances. 

The results obtained from the two-region model are in good agreemert 
with experimental observations. In particular, it is shown that the effect 
of the proton thermal anisotropy upon the angular motion is small and does 
n^ t significantly increase the predicted values for the azimuthal velocity 
at 1 A.U. 

From the solutions obtained for the models described above, we find 
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that the amount of magnetic field energy converted into kinetic energy 
in the solar wind is small and has little effect upon the expansion 
velocity. The most important effects of the magnetic field are the 
retarding torque exerted upon the outer layers of the sun’s atmosphere 
and the inhibition of the flow of thermal energy across the field lines, 
making possible a more complete conversion of this type of energy into 
kinetic energy in the acceleration region. 

Finally, we find that in order to obtain reasonable values for the 
p^article density, magnetic field intensity and energy fluxes at 1 A.U. , 
it is necessary to use a reduced value of the thermal conductivity co« 
efficient. 

This value is approximately 1/6 of the classical Spitzer’s value 
and leads to coronal densities which are almost two orders of magnitude 
lower than observed. These results give support to recent theoretical 
work indicating that the magnetic field and plasma instabilities play an 
important role in modifying the plasma transport coefficients. 




\ 
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APPENDIX A 

The computer programs developed to obtain solutions to the one- 
fluid model equations are written in APL/360 language and listed below. 

The outward integration program is "MAIN" with subprograms ’’PARAM", 

"''START'", '''S0LWIND2" and “DER”. 

The corresponding program for the inward integration process is 
'“IMAIN"" with subprograms, “TABAMl", "*PARAM!“, "START”, "S0LWIND2" and DER . 
■^e asymptotic solution is obtained by the program "ASOL” with subprograms 
"PARAMl" and "PAEAlf'. 

The density, heat flux, magnetic field intensity and kinetic and 
magnetic field energy flows in all regions, are computed by the program 
'''DBQ" o 
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VAMr//CD]7 

V MAINilCiFGiXlfiKMiKL 

[I] 'EllTER BETASTAR* 

C 2 ] BSTR^U 

[3] KL-*-KM-^Q 

[4] *E11TER ZETA' 

[5] ZETA^U 

C6] * Ell TER PllISTR IB DEGREES* 

[7] PHI5ri?-<-0x (02^360) 

[8] * ENTER GAMMA* 

[ 9 ] GAM^U 

[10] 'ENTER a* 

[ II ] 

[12] 'ENTER STEP SIZE, MAX, Z AND PRINTOUT INTERVA 
L * 

[13] IC^D 

[14] * ENTER DGAM AND NO. OF EXECUTIONS* 

[15] AL*-U 

[16] . ItAf-l 

[17] XM*-1 

[18] Mil -.PAR AM 

[19] START 

[2 0] M12'.SOLWIND2 

[21] FC-(-( 7, (p2F) )p (ZFxPSrTeaeOOOOOOOO) »( KFxt/ST), (f/Fx 

DEL^UST) ,{.THETA.F^TST) ,DODZF , (P//JFx 360^0 

2 ),BETAF 

[22] TSOL^ZFLp , ZF] . FFC p , FP] . TRETAFl p , THETAF'\ 

[2 3] MAT-^ 2 1 >S)FG 

[24] -*-(AM = 2)/0 

[2 5] -^(AM = 1)/Afl3 

[26] KM^l 

[27] ->A/14 

[28] W13:FZ;->-l 

[29] W14:-*-((FF = l)x(FAf=l))///lS 

[ 30 ] "*M 1 6 

[31] W15:/L[1]^/1F[1]«2 

[32] W16 :->-(MA = l) /M17 

[33] GAM^GAM-ALL12 

[34] -*M1B 

[35] M17:GAM^GAM+AU11 

[36] W18 :IJ-^IJ’+1 

[37] -*{IJ=ALl2l)/Ml^ 

[38] ^Wll 

[33] M19 ’.* PROGRAM EXECUTED TINES -.WISH TO 

SAVE?* 

[40] -►Wll 
V 


\ 
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VPAFAMlDl'7 

V PABAMiAiB;0ME',Gl1SfJiXPiBP 
[1] 0ME-<-2 .92E~& 

I2l GMSN^l .Z2E26 

[3] xi^i(i2)-i-iBSTR) + iUi^2)-HBSTF)*2)-{2yi20PHISTE)* 
2)iBSTR)*0 . 5 

[4] UST-^{GMSl}^0MEiGAM^ZETA)*^3 

[ 5 ] FST^ii GMSN>^ZETA *2)i GAM><0ME*2 ) * t 3 

[6] PSI^ZOPHISTF 

[7] DEL-t-PSI+ZETA 

[8] MU^{l-XI)Hl~(.XIi(.i2oPFISTF)*2))) 

[9 ] 5IG-^( l + (PSI'>il-MU)iZETA)*0 . 5 

[ 10 ] A^2 + iXI>^DEL*2)'i-(2*DEL^XI^MUxPSlil-MU)-GAM*XT 

C 11] ( 0 . Sxl+DEL*2 ) + iHU'xPSI*2 ) + { 5 ^ 2 xjTI ) -H ■i-GAMi'DEL 

yiMU^PSI 

Cl2] ALPH->-A^B 

Cl 3]. Wi455-f-(18 36x9 . 1 066£’"2 8 ) + ( 0 , 05 x 

6 .6442F”24) 

[14] BC-<-1.38E~16 

[15] TST-^(MASSt<UST* 2) i{2xBC^XI) 

[16] XP^O ,(~2xH) A2xMUyZETA*2) 

[17] VINF-^-iCUBIC ;fP)[l] 

[18] 'GAMMA: FSTAF: ' \{FST* 

69600000000);’ BETAS7AF: ' %BSTF 

[19] ’ PHISTAF: PBJSTPx ( 360^02 ); ’ XI: ' \XI 

[20] * H: ’;P;’ VSTAF: ' \UST 

[21] 'ZETA: ' iZETA\* PSI: ' ;P5I 

[22] 'DELTA: 'iDELi' MU: ' iMV 

[23] 'SIGMA: ' -,SIG 

[24] ' ALPHA: 'iALPH;' TSTAF: ';TSTi' V 

INF: ' UVINF^UST) 


V 
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V5ryli?2’C0]V 

7 STARTiA\\A2\Bl ;B2 ; Cl ; C2 ; C3 ; C4 ; C5 ; Z?1 ; D2 ; P; ffl ;Q ;F1 
\F2 iR\S \T 

[ I ] VOdZST-*-2-^{XI*DEL*2)^i2y.DEL^XIxfilJ->^PSIil-MU)-GAM^ 
XI 

[ 2 ] / !•«-{ ( 3 X WC ) - 1 + ( 1 +WC )^SIG* 2)i(l -MU ) x ( “ 1 +SIC* 

2 ) 

[3] A2<-H1-MU) 

[4] Bl-HZPr^+/llxP5J • 

[5] B2*PSI^1+A2 

[6] Clx-l+WCxPSJ*2 

[7] C2<-(DEL*2)-DEL-<MU*PSI 

[8] CZ^{2*MU^PSI)-DEL^MU 

[9] C4-(-542x;fJ 

[10] C54-C 2xPPLx//yxPSJ)+Cj4//-2xWCxPSI*2 

[II] Plt-(XLPffx;fJxCl ) + (52xC2x/lLP//x^l4PPZ,)+i42x (?-.-( >tLPPx 

XIx^CZy-PSI)^DOOZST^2>-{loP}!lSTR)*2 ) 

[12] P2x-(i4LP//x;fIxC2xBl4DEj[;) + (yllxP) + ( ( {ALPR^XI^C^)- 

2 )^D0BZST)^{ALPH^XI'^Cl)-2 . 5xPCPZS?*2 

[13] Plf-( (C'*- 2 xDPf.xZJx/fCxPSlT( l-/fC)*2 )x(( 1 -//£/) x (>12+52 4 
555) )-WC)+(2xDFLxJIxS2)-51 

[ 14 ] 52-^( Qx ( f 1 -MU ) X ( >3 1 + 514 555 ) ) - 2 ) +50DZ5T+ ( GAM^XI ) + ( 

2x5P£xJirix5l)-52 

[15] Fl^i2yXI)-(R^( XlyMU^PSI* 2)4(1 -MU ) * 2 ) x ( 2 x>l 2 - ( /?Cx 
X2)+Wy)+’l 

[16] 52^-5C5ZS2’+2x5x(>llxl-W)-l 

[17] 5-<-(51-F2 ) 42x51 

[18] r-<-((S*2)+52451)*0.5 

[19] 575Z5T1-H5tr 

[20] 5F5ZSZ’2'«-S-2’ 

[21] '575Z5T1 * iDVDZSTl 

[22] ’5F5Z5!T2 * ‘,DVDZST2 

[23] *DODZST ' iDOPZST 

[24] FC^l+PFPZSrixO .0001 
[253 TPPT/IS-^l+PCPZCrxO .0001 

[26] ZC^-l + O.OOOl 

[27] XS^iZS 

[28] 3’P^Ci4W,B5?'5,Z52’/l,5,P5-r52’P 

[29] TfI^MU,XI ,DEL,SIG,ALPH^UST,TST,RST 


7 
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VS0LWItW2iUl'^ 

V S0LUnW2 ;.Y;y;X;P;I;.7';F71;Zi';Hl;i?2;J?3;^4;/?5;..6 
L] IVl^KS , VS , THETAS 

) ] zF-^VF-<-lIF^TnETAF^VODZF->-?niF’>-BBTAF^\0 
3] Y-ej;^l[l],ICCl3,(-ICt2:) 

+ ] y4-m[3] .ii'iC 2] 

5 ] I^ZS 

5 ] Rl->rALPH'^KI*(.20PHISTR)*2 

7] R2^XI^ZETA*2 

8] RZ^ZETA*2 

9 ] Rh<-GAM^XI 

10] C25P:Y-^(4 , p , M-*-Y 

11] -*-((L[l] = 0)><(L[2] = 0))/(7J-p 

12] ■*{LllZ<0) /CEP 

13] £7£’P:-»-( J<3 + pppY[e; ;]-<-Jt-^( A'C 1]+P) PEP. W^7+Y[ 

"1 + J-^«^+1;]*P-»-XC2]t 2 2 1[«7])/C’EP 

14] PK1^(Y[1]^+/Y[ t2]),y-y+(;fC2]^6)x + /Cl] KL 
122334;] 

15] CKP'.-*-{iiFVill'i)^I) /CFP 

16] COP:-^-((Jr[2]xY[3]-YCl])>0)/CPP 

17] - MA^2 

18] -*CGP 

19] CFPiZF^ZFA^FVlill) 

20] I-<-ZP[p.Zf]xl + ICC3] 

21] 7F^-l'F,fVlC3] 

2 2] THETAF-^THETAF,FV1L21 

23] DOPZF-^PODZF ,(-El^'^^EVlllJ*2) 

24] ■*iipD0DZF)^l)/MHA 

25] -*■(( \DODZFIpPODZFD>(\DODZFL(.pDOPZF)-11))/CJP 

;2 6] MHAi-^COP 
.2 7] CRPi'PODZ PCS A 
:2 8] AM^l 
;29] -’^C’GP 

:30] CJP: ' TEMP .<0 OR DOPZ lECREASIHG' 

;31] MA^O 
"323 ~^OGP 

[33] CGPtPHIF-*-(.~2oi {ZETA^ (.SIG*2) -ZF*2) i (.VF^ZF) -MIHZF) ) 
+ ol 

;34]' WF-*-i(.VF’^3oPEIF) + ZETAy-ZF)^DEL 

135] BF!Ti4F'^( (P5TPx( l+p5J*2)xZF*2 )^l'Fx(l + ( SoPffJF)* 
2))>^TEETAF 


7 
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VD£'i?CD]v 

V Zf-jr VER \YX\Y2xT\PiAiB\C 

Cl3 -*■(//[ 1]<0 )/Dl 

[ 2 ] P^-(,ZETAiMl.2']y-X)^<.A^l-X^X*SlG*2)->iB^iX-MU>^Xx-XiMi 
2 ] 

[33 Wl*e-(/?lxl + PxP)4/f[l3*2 . 5 

[4] W2-f( 0 . 5x(W[2] - 7Ji';P)x/f[2 3+I'I//F) + ( 

2 .5x;?[l]TXI)-(Gyl.'-^x2r)-C-^0.5xP3xZxArxBxBx((pj(7* 

2)-MUHli2l)*2 

[5] M2«-W2+//yxP3x(^xB4/r[2] ) - i m'P 

[6] yi4-{2 X Cx;rj) + ( ( 2xWi/xP2f//[2] )x(ylx/txB*3)-/lxPxB) + (2:x 

WlxW2) + (2x/f[l] )^2rxP4 

[ 7 ] y 2>^ ( 2TI X w[ 2 ] x/.;[ 2 ] ) -;/[ 1 3 +P 2 *MU*A x i4 x ( P * 3 ) f W[ 

23 

[83 Z-<-(il/lx//2 ) , (-71xW[23Tyxy2) 

[93 -^0 

[103 P1:Z+- 0 0 
7 


vrBC[037 
V DBQ 

[13 'ENTER K' 

[ 2 3 

[33 NST-'riALPH'X'Ky- ( TST*2 . 5 ) x ( 20PPI52’P ) *2 ) i2xUSTxRSTxBC 

[43 NF-^NSTiVF^ZF*2 

[53 BP-^( (il/yxo4xW455x775!rx£/S2’*2 )* 

0.5)^(ZF*2 )x(2oPPJP) 

[63 QF^-K*((.THETAF^TST)i<2 , 5)x.DODZF^TSTiRST 
[73 PRF-^i UZF^RST^BF)*2)x(USTxVF^iloPHIF)*2)-DEL><UST 
X ^^px ( 1 OP//IP ) X { 2 OPPIF ) ) * 04 

[83 KEF-*-iMASSxNFxVF^UST^(ZF^RST)*2)^ 

0 .bx(iVFxUST)+WFxDBLxUST)*2 
[9 3 F^RxNSTxUSTx (,RST*2)xMASSxUST*2 
[103 'TOTAL ENERGY FLUX' iF 
[113. NBQO*- 2 1 <si(6ApZF))p(ZFxRSTi 
6 96 00000000 ) ,//F,PF.C?F»PPF,FFF 


7 
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[13 

[23 

[33 

[43 

[53 

[63 

[73 


[83 

[93 

[103 


[13 
[23 
[33 
[43 * 
[53 


[13 

[23 

[33 

[43 

[53 

[63 

[73 

[83 

[93 

[io3 

[113 

[123 

[133 

[143 

[153 

[163 

[173 


VJM>lI7/[n37 

IMAn]\FG 

'ENTER STEP SIZEy MAX, 2 ANN PRINTOUT INTERVAL' 

IC*-U 

PAR AMI 

PARAM 

START 

SOLWIND2 

FG^(,7 , (pZF) )p (ZFxT?^!?^ 696 00000000 ) , (VF^UST) ,(UF^ 
NEL*VST) » ( THETAF'>^TST) tDONZF ,{PUIF>i2&Qio 
2 ),BETAF 
MAT-^ 2 1 
MAT 

-*-(/W = 2)/0 


7Pilffi4«l[D3 V 
PAR AMI 
GAM-rTPi 13 

BSTR^TPi2^ 

ZETA-^TPiZ'\ 

R^TPW 

PUISTR-<rTPi^'] 


NPARAMIUIN 

PARAM ',A',B% OME ; GMSN iXP\BP 
0ME*-2 .^2E~& 

GW5//-»-1.3 3E2 6 

XI^{ ( 4 2 ) + 4S52’7?) + ( ( ( ( 42) + 4B5!Ti?)*2)-(2x(2oPPI5rj?)* 
2 )4BS2’P)*0.5 

UST-<r ( GMSN y^OME^GAM^ZETA)*i2i 

RST-<-{{GMSN^ZETA*2)iGAM>i0ME*2)*iZ 

PSI^ZOPUISTR 

DEL^PSI^■ZETA 

MU*-{l-XI)^{l-{XIi{i.20PHISTR)*2))) 

5J(7-^(l + (PFJxl-Wy)42F?>O*0.5 

A-<r2-y{.XI-<DEL*2)^i2^NEL-x.XI^MU*PSIil--MV'i-GAM^XI 
B^-XJx(0 . 5xl+DPL*2 ) + (WyxPSI*2) + ( 542x^1) -B+G/7/+PPB 

xMU^PSI 

ALPH^AiB 

W45S'«-( 18 36x9 .1066E~28) + (0.0 5x 
6,6442P"24) 

BCfl.38P”l6 

TST*-iMASS^UST*2)i(,2yBC^XI) 

JCP^O, (-2xB) , (2x//£yxZP!ri4*2) 

VINF*-iCUBIC ZP)[13 
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Cl] 

C2] 

C3] 

[4] 

[5] 

[ 6 ] 
C7] 
[ 8 ] 
[9] 
CIO] 
Cll] 

C12] 

C13] 

C14] 

CIS] 

C16] 

C17] 

C18] 

C19] 

C20] 

C21] 

C22] 

C23] 

C24] 


VSr<4^?2’Cn]V ^ 

START \A \ ;/l2 ;B1 ;B2 ;Cri ;C2 ;C3 ; C4 ; C5 ;P1 ;P2 1 


*TQ.*R*S*T 

hoDZST-^2-i-(.XI*DEL*2 ) + ( 2xPFix;?Ix//[/xP5l4 1-WJ/ ) -C/l/fx 


3xii/y)-l + (l+Wy)x5JC*2)T(l-Wi/)x( l+5Iff* 

2) 

A2^-i{l‘MV) 

Bl-^ZETA-^rAl^RSI . 

B2-t-P5Ixl + i42 
Cl^-l+//yxP5I*2 
C2f-{DEL*2)-DEL^MU^PSI 
C2,*-i2^lWy-PSI) -DEL^MV 
C4-t-5^2xA:i 

C5^(.2xDEL^MU^PSI)+GAM~2^MU^PSI*2 

D1-^{ALPH’<XI’(C1 ) + ( B2xC2xylLP//xXl4PPIr) +X2x ( P-t-(i4LPff’' 
jfJxCS xPSI )+DOPZS2’x2 X ( loPfflBTP) *2 ) 

D2'*-(.ALPH^ XI^C2xBl^DEL) ■^(Al^P)i‘i(. (4LPBx AfJxC4 ) - 

2 )xDODZST) + (ALPH^XI^CS) - 2 , 5-xDODZST*2 
El^(iQ-^2^DEL^XI^MU>iPSI*(l-MU)*2)y(.( 1-MU)^(A2+B2i- 
DEL))-MU) + (2>^DELy^XI->^B2)-Dl 

E2*-iQy-iil-MU)^iA 1+PliVEL ) ) - 2 ) + DOPZS2’+ ( CzAM^XI ) + ( 
2^DEL*XI^B1)-D2 ^ ^ , 

Pif-(2xxj) -(p^(;irjx/fyxp5i*2)4U-wr/)*2)x(2x>i2-(Wx 

/i2)+wy)+i 

P2t--BOBSS3’+2xPx{/llxl-//{/)-l 

B-t-(El-P2 )^2xFl 
T-*-({S*2)+E2iFl ) * 0 . 5 
DVDZSTlf-S+T 
DVDZST2^S-T 
VS*-1-DVDZST1^0 ,0001 
THETAS-^l-DODZST^O .0001 
ZS-f-1-0 . 0001 
XS->^iZS 
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^S0LWinD7WlV 

V S0LWJl]D2iX\Y-,K',D\l\J\FVl‘,LiRl \R2\E2 ;/?4 ;i?5 ;/?6 
Cl] IV1*-XS ,VS ,'HBETAS 

[ 2 ] Z F-«- VF*-BF-*- THETA F*-DOVZF-^PHIF*-BETAF-^ i 0 

[3] ^f-e-J’/lClD.JCC 1] , (^ICC2]) 

[4] y-«-iriC3] ,IF1C2] 

[5] I-<-Z5 

[ 6 ] Rl*-ALFH^XIy^(.20FHISTR)i<2 

[7] E2<-XI-<ZETA*2 

[8] R2^ZETA*2 
C9] Ri^^GAM^XI 

[10] CDP\K^{^,p,X)pL^XU^11 DER M^X 

[11] CffP 3 + PP pif[ J;3*<-<;i^U[ 1]+P) DER 

“l+Jt'J+l;]xP-fX[2]f 2 2 lUD/CEF 

[ 12 ] -»•( (L[l] = l)x(L[2]=0))/C/iP 

[13] -»-( a[l] = 0)x(L[2] = l))/CJP 

[14] P7i^(i:[i]^+/-Y[ i2]) , j-^j+(;r[2]^6)x+/[i] xi 
122334s] 

[15] •»■( (tPP'ICI] )<I)/CPP 

[16 j COP;->((J:[2]xJr[33-Z[l])>0)/CPP 
[173 MA-^2 

[18] -^CCP 

[19] CPP:ZP>hZF, ( tF71[1] ) 

[20] JvZF[p,ZF]-JC[33xZF[p,ZF3 

[21] 7F-<-l^F.F71[3] 

[22] TRETAF^THETAF ,FV1121 

[2 3] DODZF^DODZF ^ ( -Z[ 1 3 *F71[ 1 ]*2 ) 

[24] MHA:-*C0P 

[25] CffP: 'DVDZ<0 DEHOM,>0* 

[26] -^CCP 

[27] CtTP; ’P7PZ<0sFWf.>0, DEH0M.<0> 

[2 8] CGP!PffIF-<-("3o((ZF2’i4x(SIff*2)-ZF*2)T(7FxZF)-WyfZF)) 

+Ol 

[29] P/F-<-((7Fx3oPffIF) + ZFr/.xZF)TPFL 

[30] BETAF-^( ( BSTR^ ( ltP5J*2 ) xZF*2 ) ^7Fx ( 1+ ( ZOPEIF)* 
2))^THETAF 
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[ 1 ] 

C2] 

E3] 

[4] 

[53 


[ 6 ] 

[7] 

[ 8 ] 

[93 

[103 

[113 

[123 

[133 

[143 

[153 

[16] 

[17] 


[13 

[ 2 ] 

'[3] 

[43 

[5] 

[ 6 ] 
[7]. 

[83 

[93 

[ 10 ] 

[ 11 ] 


Z-^X DEE MiMliM2;YXiY2iT;PiAiBlClQ 
<2-eixi/x( J>2 ) 

->(W[13<0)/Z)1 

P-<— ( Z F ryl f //[ 2 3 X I ) X ( >1^ 1 - Zx ir X SIC* 2 ) X S'*- T 1 - Wi/ X z X X 4M 

2] 

wi-<— (PixitPxp) ?;/[ 13*2 . 5 

W2-<-( 0 . 5X (W[2 3-ri/;P)xW[23 + l'IA’P) + ( 

2 .SxMilliXI) ~(GAM^X)-C*-0 . 5xP3xIx Jxpxpx ( (SIC* 

2)-MCtW[2])*2 

Af 2 -(-M2 +/f C X p 3 X ( yl x-s ^ W[ 2 3 ) - -s- Vlff P 

Il-f-(2xCxXI) + ( { 2x//i/xP2T/./[23 )x(/!x>lxS*3)-y!xBxS) + ( Jx 

WlxA/2) + (2xW[l3 )-ATxP4 

Y2't-(IIx//[2]xA/[2] )-Af[l]+P2xWUx/lxylx (B*3)4Af[ 

23 

Z*-(A/lxAf2 ) ,. ( -YlxA/[23fIx Y2) 

-v( (yi£ 0 )x(J 2 < 0 )x(C=l))/D 2 
P3:-»-((yi^Y2)>0)/0 
Z^{Y1<0)»(Y2<0) 

■*0 

Dl'.Z-*- 0 0 

P2:'P1 • ;( tADxPSlTSOeOOOOOOOO 

C-<-0 
-*P3 


VPBC[D]V 

DBQ 

CENTER P’ 

K^O 

NST-^iALPH^K>( { 151*2 . 5 )x ( 2opHISTR)*2 ) i 2xUST^RST^BC 
NF*-NSTiVF^ZF*2 

BF*-{UiU^oH^MASS^NST'xUST*2)* 

0,5)^(ZP*2)x(2oPSIP) 

CP*-Px( {THETAF^TST)*2 . 5)^D0DZF>^TSTiRST 

PRF^i ( (ZF^RST^BF ) *2 ) x ( SSI* ppx ( loPSIP ) * 2 }-DEL^UST 

xR'Px ( l 0 PPIP)x ( 20PHIP) ) r04 

PSP'e(/.f/5SSxPPx7PxCSrx(ZPxPSr)*2 )x 

0 . 5x(( VpxySI)+f>'PxDPLxSSr)*2 

F*-H^NSTxUST^(RST*2 )x/fylSSxCSr*2 

*fOTAL EllERGY FLUX' iF 

RBQI^ 2 1 <s?(6 ,(pZP) )p(ZPxP5I^ 

69600000000), A7P, BP, OP, PRF , KEF 


V 
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v/isoi[D]v 

7 ASOLiPiSiA 

[1] PAH AMI 

[ 2 ] PAR AM 

[3] Pi-( 1 . 5x7J//F*2)-fl 

[4] 5-f-a . 5 xKI//F*2)+// 

[ 5 ] ClO^GAMiP 

[6] A-^(TSOLl2l-(.Xl^C10x(SxC10)~GAM)xTSOLll'i* 

~2 ) + (T50LC 1] *~44 3 ) + ( (rFOiiC 1] *“74 3 )x ( 7’^y.GAM)- 

35xFxC10) r9xp 

[7] C'll^-5xi442xJ!fJxP 

[8] C12^-0 

[9] C13>f--( 5xC10x(5xC’10)-G/J/)42xp 

[10] £714-*-( ( 15xG4Ayx/IxC’ll)-35xC10x(ZIx5xCll )+i4 ) 4 
9x;flxp 

[11] (721^0 

[12] C2 2-^(;SfIxC10x(FxCi0 )-<7/1A/)4;} 

[13] C2 3-t-( ( 14xC10x(;irjx5xCll)+/O-6xCylWxj:jxC’ll)4 

9x>l 

[14] F7FC-^-l, CIO, Cl l.C12, (713, <714,0, 0,0 

[15] r7FC-<-0,0,l,C21,C2 2,C2 3,0,0,0 

[16] Pr2V’FC-^0,0,0,0.0,4,(5xC21),(6xC2 2 ),( 7x(72 3) 

[17] 'ENTER MAX. Z AND PRINTOUT INTERVAL' 

[18] 

[19] Z F*-VF^WF*-TRETAF-^D0DZF*-PEIF->-BETAF*- 1 0 

[20] F3’-<-2’SOL[l]*-43 

[21] CCi:FP7/'^i ,F?*2 + X 8 

[22] V-^*/VINFy^VVEC><EPN 

[23] r-<-+//lx!rP'FCxFPi(? 

[24] Pr*-+/(-/43)xPT27FCxFP/'/ 

[25] ZF*-ZF, (Fr*’3 ) 

■[26] VF-<rVF,V 

[27] DODZF-rDODZF ,DT 

[28] 2’ffF!T/lF-«-2’PFr/lF,r 

[29] Fr^(ZF[p ,ZF]xl+J7[2] )*-43 

[30] •+(ZF[p ,ZFlsI^[l] )/CCF 

[31] PPJF-f-(“3o( (ZFT/tx(5JC*2)~ZF*2)4( KFxZF)-MC4ZF) )401 

[32] ^/F^( (7Fx3oPPJF)+ZF2’/lxZF)4PFL 

[33] BETAF^i ( BSTR x ( 1 + PS I* 2 ) x ZF* 2 ) 4 FFx ( 1 + ( 30PP JF ) * 
2))y^THETAF 

[34] MAT*- 2 1 6)(7 , (pZF) )p (ZFxP5r4 

69600000000 ) , ( VF^UST) , (l/F^DEL^UST) , iTHETAF^TST) , 
POPZF, (P//IFX 36 0 402 ) ,BETAF 


V 
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VPARAMllUl^ 

V PARAMl 

[ 1 ] GAM-^TPi 1 ] 

[2] B5ri?^-7P[2] 

[3] ZETA^-TPl^l 

[4] //«-rpC4] 

CS] PHISTR^TFiS'l 

V 

vp/ip/iw[n]v 

V PARAlUAiB-, OME ; GMSR iXPiBP 

[1] 0ME-<-2 .S2E~& 

[2] GMSN->-1.3 3E26 

[3] (f2) + ’BS?P) + ( ( ( (T2) + ^B5rP)*2)-(2x(20PPI52’P)* 

2)iBSTR)*0 . 5 

[ 4 ] UST-<-i GMStl X OME iGAM»-Z ETA ) * ^ 3 
C 5 ] RST*-{{ GMSN^ZETA *2)i GAM^0ME*2 ) * i 3 

[6] PSI-^20PHISTR 

[7] DEL^PSIfZETA 

[8] A/y-^( 1-jj) T(i-(;fi4( (2oppi5rp)*2) ) ) 

[9] ' SIG^<-(.l■^(PSI>^l-MU)iZETA)*0 .S 

CIO] A<-2 + (.XIy^DEL*2) + {2^DEL>^XI^MU^PSIil-MU)-GAM>fXI 

[11] B*-XI>^(0 .5^1 + DEL*2 ) + {MUitPSI*2)&(5i2xXl)~H+GAM+DEL 
xMUyPSI 

[12] ALPn-^AiS 

[13] W/155x-( 18 36x9 . 1066P"2 8 ) + ( 0 . 05x 
6 ,6442P"24 ) 

[14] P(7-<-l .3 8£’"16 

[15] TST->-UiASS*UST*2)r(2xBC^XI) 

[16] XP^0A-2^R) A^xMUxZETA*2) 

[17] VINF^^^iCUBIC XP)[1] 

V 

VPBQ[0]V 

V DBQ 

[1] ^EETER K' 

[ 2 ] 

[3] NST->-(.ALPI!'xK^{TST*2 .S)^(2oPRISTR)*2)r2y(VST^RST^BC 

[ 4 ] ‘ NF-i-NSTiVF^ZF*2 

[5] BF^{ (AfPxo4xW/l5'PxP5rxt/5r*2 )* 

0.5)^(ZP*2)x(2oPffJP) 

[6] QF-«--ifx ( (TPP!r/tFx2’52’)*2 . b)^DODZF^TSTiRST 

[7] PRF-*-((, (ZPxP52’xflP)*2)x(yS7xKPx(loPffIP)*2)-PPLxy5r 
X f/Px ( loPHIF ) X ( 20PRIF) ) x04 

[8] KEF-^{MASS^NF^VF^UST^(.ZF^RST)*2)^ 

0 . 5x( ( FPxy5!T)+f/FxPFLx:/5r)*2 

[9 ] F-*-H^NST'xUSTx(RST*2)xMASS’<UST*2 

[10] 'TOTAL ERERGY FLUX' -.F 

[11] UBQA^ 2 1 ISJ(6 , (pZF) )p(ZFxi?52’v 
69600000000), RF ,BE ,QR ARE AEF 


V 
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APPENDIX B 


The divergence of the anisotropic proton pressure tensor 
may be obtained as follows. In index notation we may write (Bl) as 


?2l= 2.CpaS,:y+ i. 


(B2) 




where 6.. is the Kronecker's delta. 
It follows that 




(B3) 


but from Maxwell’s equations 
hence we obtain 


(B4) 




(B5) 


which is the desired result- 







167 


APPENDIX C 


The numerical integration process is optimized by expressing the 
system of governing equations in reciprocal space. Thus we Introduce the 
new independent variable X, where 

ys-z"* 


and it follows that for any function we have 

The system of equations (III. 2. 25) then becomes 

bi. ^ •+ ’’‘'I 


(Cl) 


(C2) 


dx 


(C3) 


ax ^ ax dx 

where the coefficients b_, (j=l,8) are given by the following expressions: 


b,2. = - C^h) ] 

b.» i b|4 = VsinV/z 


(C4) 


(C5) 


(C6) 


b,6 = b.«. =o •. *>'V = v/i 


(C7) 
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t>.« - -4- 


(C8) 

(>^j — J. (C9) 

'oix -- ~ [Bv^+^CeH-ex^^os-Zj^coi^^-^Vx®] (cio) 

• bii = _ dolj^ sinj^ J {>2i(=- fc>2^ (Oil) 

^2» 

bjs = bzi = t’lT =o 

t>rt = _ Ji!_(e„-ex)sm?/<-o»V_ s| iwv^ (C13) 


bj, = 2co»^rg^^,jj^^(„30j^£,i„'5,j/4'g,^3_4a„<asV/^J 

V L (C14) 

^ <;o5^ S«rt ^f6n(l-Ko Qj . | v\) =. Qu J 

bj3 — I- 

b^ “ b'sg = •= b^7 = O (C17) 
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^ V St r»V/ <) )- Q 

(C19) 

ktZ = + Qj,(\-Qj^cci^/Z%^^)l (C20) 

b4J =: e» ; ^44 = — I (C2i) 


^45 * ^ 4 ^ “ ^47 


47 


(C22) 


^48 = 


fc> 5 j = 

b 
b 


'ffi = 


'S’) 


~ ^-3^-([0x/x) (^l+C6>^J (ji~ Q xC6^/z%\T) 

cosjz$ 4'<i?ji s«n*',^/x’J 
coSjji SiV»^ ^©M- ^^?ii CoS^/x*^ 

. b^4 = ; b^5 = 


VX 
I 


x^ 


bs6 -= bs7 = o 


bfei = ^i^(6x +<5x<^S5jf/x^) 

b 62. •= — <loS^ &tVi 5^ C -*. 


(C23) 

(C24) 

(C25) 


(C26) 


(C27) 


^58 => +Qfi^0S$S^^/+Cd5^^/x*^ (C28) 


(C29) 


V 


(C30) 
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bii = o -J bfe4=| ; 1 j6S=0 

(C31) 

; b„ =0 

(C32) 

= ^ ( ! 40>i^ 3 ( ©J. C05;// X*} 

(C33) 

^71 -s 1^72 =: ^73 = ^74 = ^75 ^ sQ 

(C34) 

b ;7 c. - CA^t^ 

(C35) 

by? = ^<>>^3 j^^C'^ + SV)+ ^(50aH-48x + 6j.^ + 

X ^ CjO'^^ S irtjzJ ^0j|‘ 


ygX' U*? 4- ('<5„ + <3x) Co5^/2X*J 

(C36) 

and 

(C37) 
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APPENDIX D 

The computer programs developed to obtain solutions to the two 
fluid, two-region model are listed below. In the inner region the pro 
grams listed in Appendix A are used to compute the solutions up to the 
transition point. 

The outward two- fluid integration program is “MAIN" with subpro- 
grams "PARW-f, "PAltAMl", "START", "S0LWIND2", "DER", "S0LWIND3” and "DERA" 
The asymptotic solutions are obtained with the program ASOLAN* and sub 

programs "PARAMl", "PARAM" and "COEFF’*. 

The heat fluxes, density, magnetic field intensity and energy flows 

are computed with the "NQTF" program. The program "RATIOS" calculates 
the proton temperature anisotropy and proton-electron temperature ratio. 
"DFP" computes the parameter values that determine the proton velocity 
distribution function. 

The constant contour curves of the distribution function may be 
obtained with the program "FUN” and subprograms "F", "G", "ZERO" and 
"2ER0P" . 
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7W/1Ia;[0]v 

7 MAni-,IC\FG\XM\KM%KL 

[1] yENTER BET AST AR' 

[2] B5rJ7-*-D 

[3] KL^KM-^0 

[4] ^ ENTER ZETA' 

[5] ZETA*-U 

[6] *EUTER PHISTR IN DEGREES' 

[7] PHISTR->-U^ (O2i3&0) 

[8] 'ENTER GAMMA' 

[9] GAM*-U 

[10] 'ENTER RATIO OF GAMPAR TO CAMPER' 

[11] (7-^n 

[12] 'ENTER H' 

[13] 

[14] 'ENTER STEP SIZE, TRANSITION Z, PRINTOUT INTERV 
AL AND MAX.Z' 

[15] IC^U 

[16] 'ENTER DGAM AND NO. OF EXECUTIONS' 

[17] AL*-U 

[18] IJ^l 
[193 XM^l 

[2 0] Mil '.PAR AM 

[21] START 

[2 2] M17'.SOLNIND2 

[23] FG-^il , (pZF) )p (ZFxJ?5!T-f69600000000 ) . (VFxVST) , (F/Fx 
DEL^UST) . ( THETAF^TST) ,DODZF, ( PffIFx 36 04o 
2),BETAF 

[24] TSOL-^-ZFlp ,ZF1 . FF[p , VFl,THETAFip ,THETAFl ,F/F[p ,F/F] 
,PHIFLp ,PHIF1,D0 DZFIp ,D0DZF2 

[2 5] MAT'^ 2 1 ^FG 

[26] •*(MA = 2) /M20 

[27] i^21 :->-(W>l = l )/W13 

[28] FW-f-l 

[29] ^W14 

[30] W13sFI'f-l 

[31] W14:-^((FI = 1)x(a:/^=1))//'/15 

[32] -*M1 6 

[33] W15s/lL[l]-<-AL[l3^2 

[34] M1B:-^(.MA = 1) /MU 

[35] GAM'^GAM-ALilli 

[36] -*-M18 

[37] MU '.GAM-^-GAM+ALZll 

[38] /fl8 ;I»7-^IJ'+1 

[39] -»-(r«7=i4L[2])/W19 

[40] -‘•Mil 

[ 
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[41] 

[42] 

[43] 

[44] 


[45] 

[46] 

[47] 


M19:' PROGRAM EXECUTED TIMES\VISU TO 

SAVE? ’ 


■*M 11 

li2C>‘.SQLWIllDZ ^ ^ 

F(7/l-^(ll.p2F)p (ZFx/?52’469 6 00 0 0 00 0 0 ) .(.VF^UST), 
DEL>^USTy^ZF) ATPL^TST) , (TPR^TST) ,QPL , QPR ATE>^TST) 
.FTFFZ , (PFiFx 36 0^-02 ) ,FFF/.F 

'T’SO LA<-ZFl p ZF ] , ’/F[ 0 KF] , [ 0 TPL ] . ^PFC p FPF ] . QPU P 

CPF ] , QPP [ P Q PP ] . t P rP ] ► C P 3 P ^ 

■*(.MA = 2) /O 

■*M 21 

V 


[ 1 ] 

[ 2 ] 

[3] 

[4] 

[ 5 ] . 

[ 6 ] 

[7] 

[ 8 ] 

[9] 

[ 10 ] 
[ 11 ] 

[ 12 ] 

[13] 

[14] 

[15] 

[16] 

[17] 

[18] 

[19] 

[ 20 ] 
[ 21 ] 
[ 22 ] 

[23] 

[24] 


V 


VPARAMZU^R 

PAR AM lAiBi OME ; GMSll %XP-,BP 

OME-^2 .9 2F"6 

CAfSF-^1 . 3 3F26 . 

II^(( v 2 )+-BSFP)+(( ( (’2 )+^BSrP)*2)-(2x(20PFI5TP)* 

2 )^BFfP)*0 . 5 

UST-*-{ GMSll>^0MEiGAM-x-ZETA)*^2> 

RST*-AGMSE-<ZETA*2)iGAM^OME*2)i^iZ 

PSI^ZoPUISTR 

DEL^PSIfZETA 

( 1 - .YI) 4 ( 1 - ( XI 4 ( ( 2 OPBIF TP ) * 2 ) ) ) 

SIG-^{\ + i.PSI*l-MU)iZETA)*^ A 

/5'^2+ (XJxDFF* 2 ) + ( 2xB&Fxa Jx//yxpSjTl-i''/B 'i-GAM-xXI 

B^XJx(0.5xl + BFF*2) + (WyxP5I*2) + (5J2xXI)-B + C^'^+BFF 


xMU^PSI 

ALPH^A^B 

1936x9 . 1066F 20) + (O.O5x 
6. 6442 F" 24) 

BC'^1 .38F"16 

TST^iMASS^UST*2 )i(2^BC^XI) 

XP*-0 , ( - 2xF) , (,2^MU-xZETA*2) 

VIIJF-^icUBIC XP)[1] 

• GAMMA : ' ; GAM ; ‘ RSTAR t ’ ; iRST^ 

69600000000);' BETASTARi ' ;BFIB 

' PHISTAF: • ;PBIBIPx ( 360fo2 ) ; ’ XI: 

’ H: ':B;' USTAR: ' '.UST 

^ZETA: ' \ZETA-A P5I: ' \PSI 

'DELTA: ' \DEL'A MU: ' iMU 

'SIGMA: ' iSIG 

' ALPHA: '-.ALPHA TSTAR: ';TSTi' 

IHF: 'iiVIHF^UST) 


iXI 


U 


% 
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V 


[ 1 ] 

C2] 

[3] 

C4] 

[5] 

[ 6 ] 

[7] 

[ 8 ] 
[9] 
CIO] 
CII] 

C12] 

[13] ^ 

[14] 

[15] 

[16] 

[17] 

[18] 

[19] 

[ 20 ] 
[ 21 ] 
[ 22 ] 

[23] 

[24] 

[25] 

[26] 
[27] 
[ 20 ], 
[29] 


75r/97?r[n]7 

STAFT;A1 ;/12 ;S1 ;52 ;C1 ;C2 ;C3 iC4 ;C5 ;Z)1 ;D2 ;P;E1 

;F2 ;/? ;5,-r 

D0DZST<-2+iXI-<DEL*2) + { 2xDEL^XI^MU^PSlTl-t1U ) - 
XI 

/ll-<-( (3x/W)-l+(l+«i/)x5I(?*2)r(l-Wy)x( 1+FJG* 

2) 

A2*-ii.l-JW) 

Sl^ZFr/+/:lxPFJ 
B2^-P5Ixl+i42 
{71-<-l+«yxP5J*2 
C2^iDEL*2) -DEL-x-MUy-PSI 
C3->-(2^MU^PSI)-DEL^MU 

C4*t-5f 2xAfr 

C S'*- ( 2 X DELxMU ^PSI) + GAM-2 ^MU^PSI* 2 

Dl^(ALPH^XIxCl)i-(.B2^C2xALPH^XIiDEL)+A2i(iP^(.ALPH>< 
XI-<C3^PSI)+DODZST^2x(lOPHISTR)*2) 
D2^(.ALPH^XIyC2y<BUDEL) + Ul^P)-^( i(ALPH^XI^CH)- 
2)y(DODZST) + (.ALPH>^XI*CS)-2 .5xD0D7.ST*2 
Fl-f-( (5^2xPPLx;^IxWi/xp5Ji(l-OT)*2 )x( ( l-//i/)x (yl2+P2f 
DEL))-MU) + {2i^DEE^XIxB2)-Dl 

E2->-{Q^({l~MU)>^(.Al+BliDEL))-2)fDODZST+iGAM^XI) + (. 
2^DEL^XI^B1)-D2 

Pl->- ( 2 X jsr J) - jx/./yxP5J*2 ) T ( l-//y ) *2 ) X ( 2xy}2 - ( Wx 

A2)+MU)+1 

F2^~D0DZST+2xH^iAl^l-MU)-l 

5t-(pi-F2 )i-2xPl 

2’-^((F*2)'+E24Fl)*0.5 

DVDZSTl^S+T 

DVDZST2*-S-T 

*DVDZST1 • iDVDZSTl 

'DVDZST2 ' ‘.DVDZST2 

•DODZST * iDODZST 

VS-^l+DVDZSTl^O ,0001 

THETAS<-1+DODZST^O ,0001 

ZS-«-l + 0 .0001 

XS<-iZS 

TP^GAM ^BSTR.ZETA^H^ Pill STB 

TN<-MU ,XI»DEL^ SIG , ALPH ^UST,TST,RST 


\Q%F1 

GAM* 
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'jsoLyniD2imv 

V S0LWIND2 iX;YiKiDiI;J‘.FVl-.LiRliR2iR3iR^lHSiF6 

[I] IVl^XS, VS .THETAS 

[ 2 ] ZF*- VF-^WF-^T}iETAF*-DODZF*-PHIF<-BETAF*- 1 0 
[3] (vIC[23) 

C4] y^jyics] ,jyi[2] 

[5] I^ZS 

C6] Rl*-ALRE>^XI>^i20PHISTR)*2 

[7] R2^XI^ZETA*2 

[ 8 ] RZ*-ZETA*2 
C9] RH^GAM^XI 

[10] CZJP 4 » p , 7 ) p ] DER N-*-Y 

[II] -^((L[l]=d)x(L[2] = 0))/C.rp 

[12] -^{L[l]<0)/i7PP 

[13] CEP:^( «/■< 3 + p p p^[«7 ; ] ♦-€7L'^-( X[ 1] +P ) DER /f-^Y+X[ 

"l+t7’-f-J'+l;]xp-.-7[2]^ 2 2 1[J'J)/CPP 

[14] Pl^l-«-(7[l]^+/7[ i2]),7-7+(7[2]46)x + /[1] Kl 
122334;] 

[15] * ( tPF 1[ 1] )al )/CPP 

[16] C0P:->( (7[2]xX[3]-7[1] )>0)/CPP 

[17] MA^2 

[18] -^CGR 

[19] CFP:ZF*-ZF ,{ rPVlLll) 

[20] J^2P[p»ZF]xl+IC[3] 

[21] 7P^FF,Pl^l[ 3] 

[22] THETAF^THETAF .FVlt2li 

[23] DODZF^DODZF, ( -L[ 1 ] xF71 [ 1 ] *2 ) 

[24] '<-((pP(?PZF) = l)/WF/l 

[ 25 ] -*•(( \DODZFZpDODZFl)>i \ DODZFl(pDODZF) -12 ) ) /CJP 

[26] MHAi->-C0P 

[27] CHPi 'DODZ POS . ' 

[2 8] MA'^l 

[29] -*CGP 

[30] CJP: ’FF//P,<0 OR DODZ IllCREASING' 

[31] W^-^0 
[3 2] -^CGP 

[33] CCP:PHJF-»-(~3o( ( ZFTylx (FIG*2 ) - ZF*2 ) ^- (VPxZF ) - Wi/vZF ) ) 

+ol 

[34] WF^ii VF'^ZoPEIF)^ZETA^ZF) iVEL 

[3 5] PFr^F-t-((B5FPx( ltPST*2)xZF*2)4l'Fx(l+(30PFJF)* 

2 ) ) xr//FFi4F 
V 


V 
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Cl] 

[ 2 ] 

C3] 

C4] 


C5] 

C6] 

C7] 


Z-<rX DER MiMliM7iYliY2-,TiPiAiBiC 
P<--(,ZETArMi 2l^X)>^(A-^l-X^X^SIG*2 ) 1 


2 ] 

M1-*—(R1^^ + P*P) 1]*2 . 5 

M2H0. 5x(WC2]-FJi';f )xWC2] + FJWF) + ( 

2.5xW[l]4.YJ)-(C/l//x;f)-C^0.5xi?3x;:x;fxBx5x((SIG* 

2)-WyvW[2] )*2 

A/24-M2tWyx/?3x(/lxBiWC2] )-v7I//F „ ,v 

yi-c ( 2 X CX.YI ) + { ( 2 xWi/xi?2 2 ] ) X ( /lx>l XB * 3 ) -4xBxB ) + ( ATx 


W1xM2) + {2xWC 1] )-A-xB4 

y2-*-C?TxWC2]xAf[2] )-AA[l]+B2xWyxi4x/ix(B*3)->WC 


[8] Z^(AflxW2)»(-JlxAA[2]vyxy2) 

C9] *>0 

CIO] DliZ^ 0 0 


V 
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'^SOLyiflDZlUV'J 

V S0LWIND3i XiY iAN iKiliLiJLiFViD 

[1] GPER-^(TS0Ll2'i^TS0Ll 1 ] *2 ) x ( J*0 . 5 ) x7»5f5L[ 3 ] x 2 x 

0 ,03914 713205x(2o?50L[5] ) xTSOLZ&l iALPH^ (. { 2oPIiISTR)* 

2 )x1+Gt2 

[2] GPAR-^G^GPER 

C 3 ] ' GPAR ' ; GPAR ; ’ GFER ' ; GPER 

[4] QPLl^GPAR^TSOEl3l^i(.TSOU3:i^XI)*O.S)i2xTSOU22xTSOUll* 

2 

[5] QPRl-*-QPLl>^2xGPERiGPAR 

[6 ] ZF*-VF<-WX^TPL-*-TPR-*-QPL-*-QPR*-TE-*-DTEDZ-*-BETAF‘*-PUIF-*-l)0DZF-*r7EETAF‘ 

P/F-t- 1 0 

[7] x^{tTB 0LH2) ,icii2 

[ 8 ] Y^TSOLi 2 ] , ( TSOLt'^ 2^TS0Ll 1 ] ) . {2pTS0Li 32) ,QPLl ,QPR^ ^ TSOLi 

32 

[ 9 ] I*-TS0Lil2 

[10] zILP//£’+-/!LPf7xl.04074206 

[11] AE-^{ALPHEx{20PHISTR)*2) , ( 0 , 3xXlxDEL*2) ^ (.H^XI) , (GAMxXl) , ( 
ZETAxXIxMU) » ( XIxMU) , {XIxDEL*2 ) ,(.2 xDEL) 

[12] CDP:K->-(.^ »p »Y)pL-^XlJ^l2 VERA M^Y 

[13] -^(( + /L) = 0)/C’JP 

[14] ->(L[7]<0)/C’//P 

[15] CEP',^{J<3i-pppV.\.J i2^JL-<-(.X{.l2+D) VERA M-*-Y+Ki 1 + ./-^J+l ; ] xO-«-Jsf[ 

2]^ 2 2 1U2)/CEP 

[16] FV'-^-(A'[i]^+/z[ i2]),y^y+{;sr[2]^6)x+/[ij fci 2233 
4 ;] 

[17] CKP:^(iiFVll2)i^I) /CFP 

[IB] COP :-^UX^22xX^ 32 -Xi l2)>0) /GDP 

[19] MA-<-2 

[20] ^CGP 

[21] CFP-.ZF*-ZFAiFVil2) 

[22] VF^VF ^FV122 ' 

[23] J-e2F[p ,ZP]xl + IC-[ 3] 

[24] WX^WX,FV132 

[25] 2’PL-*-TPL,FF[4] 

[26] TPR-^TPR,FVZ52 

[27] QPL^QPL,FV162 

[28] QPR->-QPR.FVi-J2 

[29] TE^TE,FVL&2 

[30] DTEDZ<<-DTEDZ . ( -Ll 7]xFK[l]*2) 

[31] ■*((pDTEDZ) = l) /MHA 

[32] -*■( ( \ DTEDZLpDTEDZ2)>{ \ DTEPZl~ l + pDTEDZ2 ))/CJP 

[33] MRAi^COP 

[34] CHP:*DTEDZ POS . • 

[ 35 ] MA^l 

[36] ^CGP 

[37] CJP: ’ TEMP<0 OR DTEDZ IRCR.' 

[38] MA^O 

[39] ->-CGP 

[4 0] (7(7P:PffIF't-("3o(ZF4VF)x(FFtxJ./2r)-ZFryl)+ol 

[41] BETAF^O . 5x((FS2’Px(l+pSP*2)x2F*2)^yFx(l+(30P;/IF)*2 ))xTF+( ( 
2x2’PP)+3’Pi)^3) 


V 


7 
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V 


[ 1 ] 

[ 2 ] 

L3] 

[4] 

[5] 

[ 6 ] 

C7] 

[ 8 ] 

C9] 

CIO] 

[ 11 ] 


[ 12 ] 

[133 

[14] 

[15] 

[16] 

[17] 

[18] 

[19] 

[ 20 ] 
[ 21 ] 

[ 22 ] 

[23] 

[24] 

[25] 

[26] 

[27^1 

[28] 

[29] 

[30] 

[31] 

[32] 

[33] 

[34] 
[ 35 .‘ 
[36 J 

[37] 

[38] 

[39] 

[40] 


\MX1 

C2FjIl-2isFJ2^1-CFJ2-<*M + TFI2-^2’FIxZ'FI^(TM[l]x.Y)x(P£'Lx/f[ 

7^)-ZETA 

52FI't-2x(SFI'^FFl2*0, 5 ) x { CFI*^-CFJ2 *0 . 5) 

Wl-«-(>li'/[l]xl + rFI2)vAf[7]*2.5 

/^2-f- ( 0 , 5 y^A-^-XI x/^[l]*2) + (>lF[2]x(/f[2]*2) • 2 ) + ( 

0.25x(5xW[7]) + (4xWC4])+W[3]) + 0,5xCFr2xWU]-A^[4] 

M2'!-W2 + (( ( C’'^PFLxS2FI^/f[l]x.Y)xW[ 2] xW[3]-//[4] )t4 ) + ((..[ 

5]+/!"[6])x |CFI^2xB)-/;;[3] + (/!i'7C4]x,Y)+/lF[5]x?FIxJ 

yl 1 l'«^7S - ( ( F-<-7l//[ 6 ] x77[ 1 ] x5 ) X J-FI2 ) +0 . 5 X (M[ 4 ] X5FI2 X 1 + 2 X CFJ2 ) + (//[ 

3] xCFl2xl-2xSFI2)' + M[7] 

X12^-(FFL + ;^)x(0.5x(/./C3]-W[4])x52FIxCFJ2)-Fx2'FI 

>ll3-»-M[l]xFFI2 + 2 

/ll4^W[l]xFFJ2+2 . / 

/17^-(7f[11^70 x(Zx7nxW2 + 2)+A/[7]+W[4] + ( 

0 0 5 xSFJ2x ( 1 + 2 xCFI 2 ) x/7[ 3] - Af[4 ] ) + (ylA/[ 7 ] x (W[ 2 ] + .Y) *2 )->!//[ 

4] x:Y 

WA:1^-411./112 ,/113,yU4,0,0 

A21^(.0. 5xF2FIxCFI2x/f[3]-A/[4])-FxTFJ 

/12 2-^-(FFL + 7r)x/l + (0 . 5 xC2FIxCFJ2x//[3]-/^[4] )~D 

/124-*-“lXi423->--A/[l]x52FI + 4 r ^ „ 

zl2 7+-((,y[nx52FIxCFI2x/fC3]-//[4])+2x;r)+/lA;[8]x/.xAf[2]4B 

//y2-<-/.21,,122,,123,;}24.0,0 

4 31^{2xCFI + //[l])x(A7[3]xCFlxl-3xW[4]x5FI244x>l)-FFl2xE4. 

4xW[5]+B ^ 

X32’^(C: + 2 )x(W[3]xl + 6x/A[4]xCFJ2 + 4x/l)-2xffxC’FI 
yl33'el-6x^/[3]xC'FI244x/l \ 

437^.( (2x.M[3]+JY)x((3x/^[4]xCFI2xl + CFI2) + 4x>l)-5FI2) + 

2xFxCFIxSFI2+Z 

A/JY3-*-/?31*/i32 ,;433, 0 , 0, 0 

/141-t-( (/7[6] xCFJxl + CFI2) +W[l]xB) + (Ar[4]x5FJ2x (AA[4]xCFJ2 + 

2x>n-l)+A7[l] 

442^(C’+2)x(/A[6]xCFr + B)+/7[4]xl-A/[4]xCFJ2+2x>l 

/144-*-(A/[ 3]xCFJ2 + 2x/t ) -1 

j447^-( (J‘A[6lxCFIx(7FI2-3) +2rxp) + (AA[4]x (1 + CFI2 )xl-W[4]xCFJ2 + 

2 Xj4 ) tA" 

WA'4-i-/l4l,il4 2 , 0 »i144 , 0 , 0 

A51^(2xCFI4W[l])x(AA[3]xCFI)+A7[5]xSFJ2+B 

A5 2^(F+2 )x,'/[3]-C’FJxF+2 

i453-^l 

/155«-2xCFI + B 

>15 7-^(2 + A:)x(W[ 3]x5FI2 ) + (A/[5]xCFJxl + CFI2)+B 
5 1 » /1 5 2 , 5 3 , 0 , /1 5 5 , 0 
A6 l-«-(5Fr2+/A[ 1] ) xF^AA[ 4]+/f[6 ]xC'FI + B 
462^-- CxF + 2 
/164-<-l 
i466-<-FFI + S 
467^( 1+CFJ2 )xF+X 
MA:6-^4 61 s,/16 2 , 0 ,A16 4 , 0,/166 

1 7 , /12 7 , /1 37 , ^4 7 , 457 , 7 
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[41] ,b )pMXl,rX2 ,MX3,nX^ J^XS,ffX& 

[42] z-<-(SL+.xj'?z7),(;nx;f2) 

[43] -fO 

[44] Z)ll:Z-<-7pO 

V 


V 


[ 1 ] 

[ 2 ] 

[3] 

[4] 

[5] 

[ 6 ] 

[7] 

[ 8 ] 

[9] * 

[10] 
[ 11 ] 
[ 12 ] 

[13] 

[14] 

[15] 

[16] 

[17] 

[18] 

[19] 

[ 20 ] 


V 


NQTF KiZllVl iPIilliZiViPHI iNSTiNFiBFiQl iQE iF iA iQL iQR iPBF iW \ 
KEF 

Zl<-MATl ; 1] x6 9600000000 f7?5!T 
BNQ^ l 0 

Vl-^MATl ;2]t{/ 52’ 
p;/llt-02x/f/ir[ ;6] ^360 
Z-t-Zl, 2P 
V<-V1,VF 

W^iUTl ; 3] ,FGi4[ 3 ; ] 

PHJ^PllIl^PHJF 

BST<r{ALPH>^K^ (TST*2 . {20PHISTB) *2) ^2^UST^BST^BC 
NF-^BSTiVyZ*2 

BF-*-( (WyxO4x/MSFx7;srx[/5r*2)*0 . 5 )t(Z* 2 )x 1 2 0P//J 
Qlir-Kx(r!ATi ;4]*2.5)x/f/.?[ ; 5]xrS2’4-PS2’ 

QE<--K^ iCTEx?3T)*2 . 5 ) x xrS2’-JF52’x 1 . 04 

F^HxBSTxUST>( (P£'^T*2 ) xMASS^UST*2 

'TOTAL ENERGY FLUX ';F;* F/R *2 *;F^( 

214x69600000000 )*2 

QL-^-QPL^Af-NST'^BC^TST^UST 

QPf-QPPx/1 

PRF^{ ( (ZxP5?xSF)*2 )x(y5!rx7x(loPBI)*2)-f/x{10P//J)x( 

20PUI) ) t04 

;fFF'^-(W/lSBx/'/FxFxy5rx(ZxP5r)*2)x0. 5x( ( Fxt/5!T)+?/)*2 

BNQ^ 2 1 fi?(8,pZ)p(2x/?S2'T69600000000),WF,BF,Cl,(?F,P,',v:,(F^( 

pQ1)po),qp,ppf,a:ff 


VP4FI05[D]F 
' V RATIOS 

[1] PF-<-D^ 2 1 6?(4,prPF)pFf;/l[ 1 ;],{Z’PBt 2’PP) , ((2’PL + 2x2’PP)xrS2'v 

3) , ( ( rpjs+2xrpp)-J3xrF) 


VFFP[0]V 

V VFP\GPL\GPR 

C 1 ] GPL-^i, 2 x(iPLx( 7 FxZF* 2 )x( XI^TPL ) *0 . 5 ) iTPL 

[2] GPR-*-(QPR^(. 7FxZF*2 )x ( XI^TPL ) *0 . 5 ) vTPP 

[3] □-*-2 1 i^{b,0GPL)0FGAi\\~^,GPL,GPRATPLiTPR)ATPL^TST) 

V 
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lAS0LAUi^'\V 
V AS0LAN\PiS 

[I] PARAin 

C 2 : PARAM 

[3] P^(1.5xKJ//F*2)-ff 

[4] Ff-(l . 5xV'I//F*2)+F 

[5] E0--TS0LAil1*-i^ 

[6] TS0LAl9^*-TS0LA\_'^'\>(fS0LAl\'\ 

[7] A->-0,A-<rTSOLAn 3 4 5 6 4 6 3 15 12 3 

[8] BOL-.COEFF 

[9] E0V-^E0*Q^\Q 

CIO] V'F-f-C’lO.C’ll.O.ClS, C14,C15,C'16,C17,C’18 

[II] TFIVI, 0,C2 2, C2 3, (724,(72 5 ,^26 .(72 7,0 

[12] rPiFi-l, 0,0, 0,0,0 ,(736,0,0 

[13] TPRV->rl^ 0*0,(74 3,(74 4,0 ,(74 6 ,(74 7,0(4 8 

[14] (?PP F-i-l.O, 0 ,(75 3 , (75 4, 0 ,C56 ,(75 7 ,(758 

[15] QPPl^-<-l, 0 *0 ,(76 3 ,(764,0 .(76 6 ,(76 7 ,(76 8 

[16] PrF*>-4,0. (6xC2 2),(7x0723 ),( 8x(724 ),(9>:(72 5 ).( 10x0726 ),( 
11XC27) ,0 

[17] P/7-*-l , 0 , 0 , 077 3 , (774 , 0 , C7 6 , C7 7 , (77 8 

[18] A2-<^TSOLAm^{EO*'^)^*/TEV^EOV 
[19 3 AZ-^TSOLAi 3]4(P(?*6 )x + /2’PLFxP(7t^ 

[20] /I4^2’F(7L/l[4]T(P(7*3)x+/!r'PPI/xP0)F 

[21] AZ^TSOLAlZ'^i{EOi>lZ)^-^ /QPLV^EOV 

[22] /16^?5(?i^[6]T(PO)*12)xt/(3PPFxPO;^ 

[2 3] 47^2’5(7L/4C9 ] f (P(7*3 )xt/P/^xp(:>p' 

[24] ->-((i42-^C2])50.001x/t2)/pyP 

[25] ^.-»-0 ,/l2 ,yl3 ,j 44 *yl5 ,/6 ,^7 

[26] -^BOL 

[2 7] BUR '.'ENTER MAX. Z AND PRINTOUT INTERVAL' 

[28] 

[2 91 *.4 VECTOR ’ 

[30] ZF'>rVF*‘NF*-BETAF-rTPL-^TPR*‘QPL-^QPR’<^TE->-DTEDZ'<rPHIF<-\Z 

[31] P2’f3’5(7L/lCl]*-^3 

[32] COL\EPN'^ET*Z 

[3 3] V-^VINF^l +PPPC 4 ] X + / VVy-EPN 

[34] 2’-H/12xPPi7[5]x + /rpp'xpp/; 

[35] 2’P4P-<-/13xPP7;[ 7] x+/TPi!/xPPW 

[36] 2’PPP'<-/14xPP/i/[ 4 ] x + /?pp7xPP/7 

[37] QPAR^^A 5 x ( EPNl 2 ] * 1 5 ) x + /QPL V^ EPN 

[38] QPER<-A & X ( EPNi 2 ] *1 2 /QPRV^EPN 

[39] D?^-(424 3) xPp;(7[8]x + /PT7xPP^/ 

[4 0] ii^^->17xPP7/[4]x+ /^/yxpp/;? 

[41] zf^zfant*~z) 

[42] VF'>rVF^V 

[43] DTEDZi-DTEDZ ^DT 

[44] TE'>-TE,T 

[4 5] TPL'^TPL,TPAR 

[4 6] TPR<-TPR,TPER 
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[ 47 ] QPL^Q?L,QT>AE 
C4 8J QPn-^QPR,Q?EP. 

[ 4 9 J WF<-WF , i/ 

C50] Er^-(ZF[p»2F]xl+J7C2] )*--3 

[51] -^(ZF[p^ZF]<FyCl])/COi 

[52] PHIF*-{ 30((ZF7’/1x(5JC*2)-ZF*2)-f(FFxZF)-/fy^ZF))+0l 

[53] BF.TAF*-0 . 5x ( (SSTFx ( 1 + P5I*2 ) xZF*2 KFx ( 1+ ( 30PFIF)*2 ) ) x^F+C ( ( 
2 xTPP)+Y'PF)t3) 

[54] □<- 2 1 6jFG/l-H(ll,pZF)p(ZFxPS3’^69600000000 ),O^FxyFT), (f/FxDFL 

^UST) . iTPL^TST) , {TPRxTST) ,QPL,QPF^ {TE^TST) ^DTEDZ , (PPIFx 
360502 )»BETAF 

V 


VP/lPi4WlCD] V 

V PXP/lWl 

[1] GAM*-TPill 

[2] BSTR-^TPZ2l 

[3] ZETA-^TP^22 

[4] - ;/-<-7’P[4] 

[5] PHISTR^TPIS^ 

V 


'7P^Pi4W[D] 7 

V P>IP/IA/;/l ;B ;C«F;F//5/;;-fP;SP 

[1] 0ME*-2 .92E~6 

[2] CWF//^1 .3 3F26 

[3] 52) + 5BF?P) + {( ( ( 52 ) + 5BSTP)*2 )-(2x(2oP//rS2’P)*2 )5BSrF)* 
0.5 

[ 4 ] U5T -*- ( CMSfl-^OHEiGAM* Z ETA ) * 5 3 

, [5] RST*-U GMS}I^ZBTA*2 ) iGAMy.0ME*2)*iB 

[6] PSI<-2oPHISTR 

[7] VEL^PSI^ZETA 

[8] wy+-(i-jsri)5(i-(A'i5((2oppisrp)*2)) ) 

[9] SIG^-(l + (PFIxl-wy )5ZF!I’>1)*0 . 5 

[10] /5-^2+(XJxPFF*2) + (2xPFLxirixWxP5l4l-wy)-G/!Wx^J 

[11] P^JfJxCO . 5xl + PFi*2)+(//FxP5J*2) + (552xZl)-p+G/lM+PFix//PxP5J 

[12] ALPR^AiB 

[13] /MSS^C 1836x9 . 1066F~28 )+{0 . 05x6 .6442F“24) 

[14] BG-^1.38F16 

[15] 2’52’-<-(AM5Sx£/FT*2 ) 5 ( 2xPGx:?'J) 

[16] XP^O A-'Z^R) ,{'2y'MUxZETA*2) 

[17] VINF^iCVBIC ;rP)[ll 

V 
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Cl] 

C2] 

[3] 

C4] 

[5] 

[6] 

[7] 

[ 8 ] 

[9] 

[ 10 ] 
[ 11 ] 
[ 12 ] 

[13] 

[14] 

[15] 


[161 

[17] 

[18] 

[19] 

[ 20 ] 
[ 21 ] 

[ 22 ] 

[23] 

[24] 

[25] 

[26] 
[27] 


VCOffFFCn] V 

COEFFiM 

Cll-(--5x/1[2] T4x/f 

C22^(2xC*10x/[iv] + UIx5xC’10)-C/l/fxJI)T/l[2] 

C13^{(4x/l[4]xC10)-/lC 3] + 5x4[2]xC22)44xW 

C23^((Cllx(6x/.[4]) + (14xZIxC10x5)-6xC>lWxXj)-;lC2]xC-10)x 

11t42x4[2] 

C’14-^-( - { 5x/t[2]xC2 3) -4xvi[4]xCll)^4xW 

C24-^Cllxl + 4x,YIxCllx5^-3x/4[2] 

C15-«— 5x/l[ 2 ]xC24 ^4xW 

C36^((4x4[5]M[ Z2) ~VINFiZETA)^VlEF^ZETA 
C^Z*--C10 

C46t(?i5*2) + (((1.5x4[6]U[4]) + FIWF^2xZE?yl)xm7F4Zffr/l)-C13 
C47-e(2xC10xCll)-C’14 

C48-f-(Cll*2 ) -C15 „ 

C25>f((2x/l[4]xC’46) + ((2xl'I//F^Z27ryl)x(yl[4]xVJ/;FTZEr;l)- 
3 xyl[ 6 ] )+2xC10x(6xA:Jx5xC’13) + (3x/1[2]xC 22 )+;![ 3]-2x>lC 

4]fC10)-V5x/lC2] ^ . , 

C16-<-( ( ( 2xyjFF^ZFrv4)x/.[6]-/.[4]xFI//Fi’ZF2’/l)-(4xyl[4]xC46 ) + 

5X/1C2] xC25 ) 44xW ^ ,, ^ 

C 2 6-»-(~7x/![4]xC4 7)t(20xC'10x{/l[4]x(7JFF*2) ^ (ZF?>1*2 ) ) + ( .YI’<5x 
£714)+/1[2 ]xC2342) 

C26-»-C26 + ( Cl lx (3x/l[3]) + (20xYIx5xC13) + ( 5x/l[4]xC10) + 

8x4C2]xC22) + (4x^[2]xC13)-C14x(5xC/inx;:J)+4[4 j 

C26-^C2649x;4[2 ] ^ 

C17-«-(("5x/lC2]xC26i4) + ( yj//F*2 )x/l[4]xCl0 4ZFr/t*2)4]>/ 
C27-^("16 x/ 1[4] xC48 ) + ("2xCllx(22x/l[4]x (7 IiVF* 2 )^-(ZFT/l* 

2)) + (22xYIxSxC14) + (”6x4[4]xCll)+9x/l[2]xC2 3) 

C2 7-^C27+(4xC15x/lL4] + 3xC4/‘/x^I)-(C10x(44xJfIx5xC15) + 

1650x^[2]xC24t61) + 4x/1[2]xC14 

C2 7^C2 7-J21x/lC 2] ^ ^ _ 

C18-»-((“5x^[2]xC2 7 * 4 ) + (CllXil[ 4 ]x( VINE *2 ) i(ZETA*2 ) )-i 4 [ 
4]xC4 8)t^M 

C5 3^-C10 + 3x4[3]x^[4]T4xXJxZF2’>lxyJffFx/l[5] 

Q 5 4-4- Cll 

C56->-(3x(yl[ 3]*2)f4xYJxZFr/!xyiFFx>lC5]) + (3x(yi?/F*2)^ 
2^ZETA*2 )-C13+C10xC10+2xC53 


[28] 

[29] 

[30] 

[31] 

[32] 

[33] 

[34] 
[35j 

[36] 

[37] 

[38] 

[39] 


C 5 7 ■!— Cl 4 + 2 X C 1 1 X C 5 3 
C58-«-(Cll*2 )-C15 

C6 3-»--(2xC10) + (;4[4l*2) i7^XI’(ZEEA^VIt}F-xAZ& 3 
C64“^“2 

C6 6'»-(”2xC10x( 2xC6 3)-i-C10 ) + (”2xC13)+4[ 3] x/l[ 4 ] 44x;fJxZF2’ylx yiFF 
x^[ 6 ] 

C6 7-^(6xC10xCll-J7) + ("'2xC14)-18xCllxC63v7 

C68^(~3xC11xC64t 2 )-2xC15 

C7 3-«-( ( C10y-MU>(ZETAiVnJF) -Al^2i2xXI^ZETA)rDEL^Al7i 

C7^<-CllxMU^ZETAiDEL^VINFxAZ72 

C76->-(C13^l}U>^ZETAiDELxVINF^AZl2)*ifWTVinF)-^(AL3li2-xDEL^ZE.A 

xYJx/l[7])-C10xC73 

C7 7^(C14x/./CxZF2’/l^DFLxyiFF) + (v1[4] ^•14xCFCxZF2’Ax;fIx^L 
7] )-(6xCllxC73i7 )t8xC74xC10i-7 
C7S^{C15^MU^ZETA^DEL^VINF^AZ7 2 )-CH«<774 
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V 


[ 1 ] 

C2] 

[3] 

[4] 
C5] 
C6] 
[7] 
C8] 
C9] 
[10] 
[ 11 ] 
[ 12 ] 

[13] 

[14] 

[15] 

[16] 

[17] 

[ 18 ] 

[19] 

[ 20 ] 


V 


^//^^^[n] V 

13QTF K\Z1\V1’, FBI 1 ; Z ; 17 ; FBI -JISI \BF \BF iQl \QE \F \A -,QL',QR FRF ; i ! ; 

KEF 

Zl*-\0 

BNQ->- 1 0 

Fl^iO 

P//Jl-*-lO 

Z-^Z1,ZF 

V-<rVl , VF 

f-/-<-,F(7-4[ 3; ] 

FHI^FHI1,PBIF 

NST*-iALPB^K>(iTST*2 ,5)^i20FniSTR)*2)i2xUSTyRST*BC 
NF-^BSTiVxZ*2 

BF->-i(.MUxoHxMASSxNST^USI*2)*0.5)i(Z*2)^ | 20PBI 
Ql->-iO 

QE^-K^{(.TEy-TST)*2 . ^)^DTEDZ^TS'nRSTy^\ 

F*-B^NST^USTy~{RST*2)->f-MASS-x.UST*2 

'TOTAL ENERGY FLUX ’ ;F 

QL->-QPL^A*-NST>^BCxTS?y<UST 

QR'*-QPR'^Ai 

PRF<-{ {iZxRST^BF)*2)x(UST^Vy(lr)PHI)*2)-W^(loPBI}x( 

20PHI) ) f04 

KEF^(MASSxP!FxV>: USTx(Z xRST )*2)x0.5x((7x PST) +P7 ) *2 

SBQA^- 2 1 5?(S,pZ)p(ZxP5fiG9£OOGOOOOO) ,;/F,PF,01.P?F,P.^i. (P^ 

ipQDpO) ,QR,PRF,KEF 


VRATIOSZn'iB 
7 RATIOS 

[1] PL*-U^ 2 1 <Sj(4,prPL)pF(77l[l;a,(2’PL*2’PP)»((rPi + 2x2'PP)xr52’4 

3) „((2’PF+2xY’PP)^3x7ff) 

V 


vPFP[n]v 

,7 DFPiGPLiGPR 

[ 1 ] GPL-^-(2 xQPL X ( VFxZF*2 ) x ( XHTPL ) *0 . 5 ) iTPL 

[2] GPR-^-(QPRx(VFxZF*2}x(XIiTPL)*0 .S)iTPR 

[3] P-2 1 H)(5,pGPL)pFGAZlil,GPL,GPRATPLiTPR) ATPLxTST) 


% 


7 
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VFy//[0] V 

V FUNiHiliJlPHIiMiXIOiK 

[1] 'ENTER GPAP.^ GPER , TPAR , TPER ' 

[2] JF-«-D 

[3] W^I//C3]4 JFC4] 

[4] XIO-^0.001 ZERO ~1 0.001 

C5] H^(*-XIO*2)^l^-UNmxXIO^i(2xXIO*2)i3)-l)-2xINL2Z^XIO 
[6 3 0 . 1 X 1 9 

[7] PAJ-^(-ol^2)+Olx(418)xil8 

[8] LP*-({pK) ppPRDpO 

[9] 

[103 COl :LP[J;tT]-^0 .001 ZEROP 4 0 
[113 J-^J+1 

[12] ■^(JipPIlD/COl 

[13] r^j+c7-^i 

[14] -v{ISpA)/CDl 

[15] LP^ 2 1 (Sf((l + pA),pPPl)p(PPIx360iO2),rIP 

[16] 'SCALE FACTOR' 

[17] D-<-5^( 2x1 . 38P"16 xJW[33 4/MS5)*0 . 5 

[18] □•t-5x(4W)*0.S 
[191 LP 

V 


VFlu] V 
V Z^F X 

[1] Z+(((4t3)x GAML 1 ] ) X A*4 ) + ( Ax 2 - GAMl 1 3 )i’GAMi 1 ] + ( 2 x GAIK 

2 3 ) - ( A*2 ) X ( 2 X G/1W[ 1 ] ) +4 x(7>l,'./[ 2 ] 

7 


VC[D]7 
V Z-*-G XlAlB 

[1] Z^A[Il-(*-(/l*2)+B)xl + (I7;Cl3x/lx(~l + 2x(/l*2)43))+2xIff[ 

2 ]x(;l^-AIO + AxloP;/J[J3)x("i + B-<-Wx (Ax20PPI[ J3)*2 ) 


VZFPO[C] V 
, V Z-^TOL ZERO B‘.T 

[1] -*-0xx2’0L2 I r^^F Z-^0 .5X + /B 

[2] -*-l,B[2l(0<F)=^0<F P]'^-Z 
7 


VZFPOP[Q]V 
7 Z<-TOL ZEROP B;T 
[11 -*“0xi2’C>Z,> I r-i-G Z-^0,5x + /B 

[2] ■M.B[2i(0<r)?0<<; Bl-^-Z 

V 
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